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applied  to  preform  design  in  shell  nosing. 
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1 .  Problem  Statement 


The  objective  of  the  research  is  to  establish  a  finite-element  formula¬ 
tion  for  the  analysis  of  flow  behavior  of  rate  sensitive  materials  under 
large  plastic  deformation.  Specifically,  the  main  concern  is  with  flow  anal¬ 
ysis  of  metal-forming  processes  at  elevated  temperatures.  The  material's 
behavior  is  characterized  as  rlgld-viscoplastic.  To  accomplish  this  objective 
the  research  was  undertaken  for  the  construction  of  constitutive  relations 
and  for  establishing  basic  principles  for  the  finite-element  formulation. 

In  order  to  demonstrate  the  usefulness  of  the  theory,  the  finite-element 
method  code  was  written  and  the  problems  unique  in  the  area  of  metal -formi ng 
were  solved. 

2.  Sunrnary  of  the  Results 

A  thermodynamic  theory  bf  deformation  using  Internal  variables  was  adapt¬ 
ed  to  rigid  viscoplastic  materials,  and  the  finite  element  discretization  for 
the  analysis  of  metal-forming  problems  was  established.  Of  particular  signi¬ 
ficance  Is  an  accomplishment  of  the  coupled  analysis  of  deformation  and  heat 
transfer.  The  finite  element  formulations  for  coupled  analysis  were  applied 
to  the  problems  of  solid  cylinder  compression  and  ring  compression.  The 
temperature  distributions  In  solid  cylinder  compression  are  compared  with 
experimental  data,  and  the  computed  results  of  ring  compression  of  aluminum 
at  elevated  temperature  are  compared  with  experimental  observations  as  well 
as  with  calculated  results  not  including  the  temperature  effects.  It  is 
demonstrated  that  the  coupled  analysis  accurately  simulates  these  processes. 

The  capabilities  of  a  coupled  thermoviscoplastlc  analysis  of  non-steady- 
state  forming  processes  were  further  demonstrated  In  the  analysis  of  the  shell 
nosing  processes  at  room  temperature,  as  well  as  at  elevated  temperatures. 

In  cold  nosing  simulations,  the  nine-node  quadrilateral  elements  with 
quadratic  velocity  distribution  were  used  for  the  workpiece..  The  treatment 
of  a  moving  boundary  in  the  analysis  of  nosing  is  discussed  and  successfully 
Implemented  In  the  finite-element  program.  FEM  simulations  of  105-mm  diam. 
shells  of  AISI  1018  steel  and  aluminum  2024  were  performed  and  solutions  were 
obtained  In  terms  of  load-displacement  curves,  thickness  distribution,  elonga¬ 
tion,  and  strain  distributions.  Comparisons  with  experimental  data  show  very 
good  agreement. 


In  the  analysis  of  shell  nosing  at  elevated  temperatures,  the  strain- 
rate  effects  on  materials  properties  and  the  flow  stress  dependence  on  temper¬ 
atures  were  Included  In  the  finite-element  analysis.  A  thermodynamic  theory 
of  viscoplasticity  based, on  rational  mechanics  was  adapted  to  a  rigid-plastic 
material  Idealization.  An  Implicit  scheme  Is  used  for  the  time  Integration 
of  heat  transfer  equations,  which  Is  coupled  to  the  plasticity  equations. 

The  nine-node  quadrilateral  elements  with  quadratic  velocity  distribu¬ 
tion  were  used  for  the  workpiece,  and  four-node  quadrilateral  elements  were 
used  for  the  die  In  the  heat  transfer  analysis  and  temperature  calculations. 

The  coupled  analysis  of  heat  transfer  and  deformation  was  applied  to  the  form¬ 
ing  of  AISI  1045  steel  shells.  Correlation  between  simulation  and  experi¬ 
mental  results  are  good.  The  results  of  the  analysis  of  shell  nosing  process¬ 
es  again  confirmed  that  the  coupled  thermoviscoplastlc  analysis  by  the  finite- 
element  method  Is  Indeed  a  significant  accomplishment. 

The  research  was  continued  to  explore  the  use  of  the  capabilities  of  the 
finite-element  technique  to  outstanding  problem  areas  In  metal  forming.  In 
the  study  on  the  finite-element  analysis  of  nosing  of  aluminum  alloy  2024 
tubular  preforms  with  uniform  wall  thickness  at  room  temperature,  the  detailed 
mechanics  were  obtained  by  computer  simulation,  and  the  finite-element  method 
was  suggested  as  a  fitting  approach  to  solving  the  problem  of  preform  design 
In  metal  forming.  The  approach  assumes  that  the  strain  distributions  computed 
by  the  finite-element  method  are  valid  for  preforms  having  nonuniform  thick¬ 
nesses  Initially,  and  two  types  of  preform  shapes  were  designed  to  obtain 
uniform  wall  thickness  after  nosing.  Experimental  results  of  nosing  these 
preforms  showed  that  the  approach  to  the  problem  Is  satisfactory. 

Most  recent  achievement  In  this  area  of  preform  design  in  metal  forming 
Is  an  Introduction  of  a  new  approach  to  the  problem  as  a  unique  application  of 
the  finite  element  method.  The  concept  Involved  In  the  approach  is  to  trace 
backward  the  loading  path  In  the  actual  forming  process  from  a  given  final 
configuration  using  the  finite-element  method.  The  concept  was  tested  with 
the  problem  of  simple  compression  of  a  cylinder  and  the  method  was  applied  to 
preform  design  In  shell  nosing. 

For  the  further  development  of  the  new  approach  to  preform  design  In 
shell  nosing,  two  approximate  solutions  were  derived.  The  degree  of  approxi¬ 
mation  was  examined  by  the  accurate  numerical  solution  based  on  the  new  approach 
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by  the  finite-element  method,  for  an  example  of  uniform  wall  thickness  of  the 
nosed  shell.  It  was  demonstrated  that  both  approximate  solutions  are  useful 
In  design,  one  for  Its  simplicity  and  the  other  for  its  accuracy.  The  re¬ 
sults  of  the  study  on  preform  design  are  encouraging  in  the  challenge  to  this 
significant  and  difficult  problem  In  metal  forming. 
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A  COUPLED  ANALYSIS  OF  VISCOPLASTIC 
DEFORMATION  AND  HEAT  TRANSFER— I 


THEORETICAL  CONSIDERATIONS 
N.  Rebelo  and  S.  Kobayashi 

Department  of  Mechanical  Engineering,  University  of  California,  Berkeley,  CA  94720,  U.S.A. 

( Received  22  February  1900) 

Seminary — A  thermodynamic  theory  of  deformation  using  internal  variables  is  adapted  to 
rigid  viscoplastic  materials,  and  the  formulations  for  finite  element  discretization  are  given 
specifically  for  metal  forming  problems.  A  scheme  of  numerical  integration  of  heat-balance 
equations  that  leads  to  the  coupled  analysis  of  deformation  and  heat  transfer  is  presented. 
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rate  of  internal  energy  per  unit  mass 
surface  density  of  heat  influx 
mass  density  of  heat  influx 
stress  tensor 

infinitesimal  strain  rate  tensor 
mass  density 

rate  of  entropy  per  unit  mass 

absolute  temperature 

specific  free  energy  per  unit  mass 

temperature  gradient 

location  of  material  point 

time 

{o^(x,  /),T(x,  t),Tj(x,  I)} 

{Mx.Ortx,  f),q,(s.  t)} 
lg(x.t)Mx,t )} 
velocity  of  particle 
body  forces 
internal  parameters 
{g(x,  t)Mx,  0} 

deviatoric  stress  tensor 

effective  strain  rate 

effective  stress 

static  yield  stress 

material  parameters 

material  constants 

energy  for  self  diffusion 

Boftzman  constant 

dislocation  density 

thermal  conductivity 

energy  per  unit  length  of  dislocation 

variation  of  (  ) 

heat  flow  through  surface  of  normal  n 
temperature  gradient  along  normal  n 
radiated  heat 
emissivity  of  the  surface 
Boftzman  constant 

reference  temperature/ initial  temperature 
temperature  of  external  environment 
surface  temperature 
converted  heat 

surface  heat  transfer  coefficient 
increment  of  time 
heat  capacity  matrix 
heat  conduction  matrix 
heat  flux 

nodal  point  temperature  and  temperature  rate  vectors 

time  integration  parameter 

intermediate  values  in  time  integration  scheme 
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1.  INTRODUCTION 

Theoretical  considerations  pertaining  to  several  disciplines  are  needed  to  perform  a 
coupled  analysis  of  viscoplastic  deformation  and  heat  transfer.  Because  these  con¬ 
siderations  adapt  classical  formulations  to  the  problems  in  metal  forming  and  because 
the  necessary  formulations  are  scattered  in  origin,  we  believe  it  important  that  they  be 
presented  together  as  a  coherent  formulation. 

With  the  availability  of  inexpensive  large  computing  facilities,  there  is  a  trend  in 
the  analysis  of  large  plastic  deformation  for  simplified  material  behavior  descriptions 
to  be  slowly  abandoned  and  more  accurate  descriptions  introduced.  The  “problem’'  of 
solving  complicated  equations  is  left  to  well  developed  computer  routines,  and  the 
analysts  can  now  devote  themselves  to  evaluating  the  appropriateness  of  the  analysis. 

The  successful  application  of  the  finite  element  method  to  the  analysis  of  metal 
forming  processes  using  rigid  plastic  materials  [1-5]  leads  to  the  simple  extension  to 
rigid  viscoplastic  materials  [6].  The  present  work  aims  at  finding  the  effects  of 
temperature  variations  within  the  body  on  its  deformation  behavior.  A  ther¬ 
modynamic  theory  of  deformation  utilizing  the  concept  of  internal  parameters 
developed  by  Perzyna[7,8]  is  adapted  here  to  a  rigid  viscoplastic  material.  As  an 
introductory  treatment,  plastic  strain  and  the  dislocation  density  are  considered  as  the 
only  internal  parameters,  and  physically  based  laws  are  proposed  to  deal  with  them. 
We  also  discuss  the  problem  of  numerical  integration  of  heat  balance  equations  and 
present  a  time-step  scheme,  leading  to  a  coupled  analysis  of  deformation  and  heat 
transfer. 

In  the  companion  paper,  some  simple  applications  are  shown,  with  calculation 
procedures  and  details  on  the  treatment  of  boundary  conditions. 

2.  THERMODYNAMIC  FORMULATION  OF  THE  BEHAVIOR  OF  RIGID 
VISCOPLASTIC  MATERIALS 

2.1  Thermodynamic  variables 

We  follow  the  work  of  Perzyna[7,8],  which  presents  a  thermodynamic  theory  of  rate  sensitive 
plasticity,  following  rational  thermodynamics,  with  dissipative  effects  described  by  internal  parameters. 

A  thermodynamic  process  is  defined  as  one  which  satisfies  two  postulates:  balance  of  energy  and  growth 
of  entropy.  The  thermo  mechanical  state  of  a  body  at  a  certain  time  t  is  a  collection  of  values  of  internal 
energy,  internal  entropy,  entropy  influx,  heating,  and  power  expended  with  its  motion. 

Using  the  principles  of  conservation  of  mass,  linear  momentum,  and  angular  momer*um.  the  first  law  of 
thermodynamics  is  locally  stated  as 


+  qw  +  pr  -  pe  «  0, 

and  the  second  law  of  thermodynamics  yields  the  Du  hem  inequality 


(1) 


(2) 


Calling  specific  free  energy  per  unit  mass  <b  to  be 


equations  (I)  and  (2)  become 


+  -e-Tv. 


»»««  +  qu  +  P»'_P(*+ Tn)*0,  (3) 

_*+f„+2fit+S(Ii*0.  (4) 

P  pT 


We  assume  that  in  order  to  describe  a  thermodynamic  process  in  a  body,  we  need  to  specify  a  set  of 
three  bask  variables 

f(x,  t)  -  {<rri«.  0,  T(m,  f),  Tj(x,  0}  (5) 


at  each  point  x,  and  seek  to  determine  at  each  material  particle  a  collection  of  functions 


P(x)  ■  (f(a.  »).»(*.  0). 


(6) 
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where  v  represents  the  set  of  variables 

ir(x,  t)  *  {gfx,  t),  n( a,  t),  <fc(x,  f)}. 


(7) 


satisfying  inequality  (4). 

Heat  supply  if x,  t )  and  body  forces  F  do  not  need  to  be  specified,  as  they  can  be  determined  by  (3)  and 
by  equilibrium  equation 


Oiy  +  F,  *  p% 


(8) 


2.2  Constitutive  equations 

For  the  definition  of  the  thermo-mechanical  state,  two  propositions  are  accepted;  (I)  the  thermo- 
mechanical  state  can  be  specified  by  g(x,  t )  and  by  its  method  of  preparation,  (2)  the  method  of  preparing 
the  actual  stress  temperature  configuration  can  be  determined  by  one  group  of  internal  parameters  w(x,  /), 
which,  in  a  ridid,  rate  sensitive  plastic  material,  describes  internal  structure  changes  produced  by  plastic 
deformation. 

Accordingly,  the  thermo-mechanical  state  of  a  particle  at  x  at  time  t  is  described  by  the  value  at  t  of 

*(*.  t )  -  <g(x,  t ),  tv(x,  f)},  (9) 

and  by  the  initial  value  problems  given  by 


tv(x,  1)  -  l V(h(x.  t )), 

(10) 

W(x,  to)  -  wrfx). 

(U) 

It  is  assumed  that  constitutive  relations  determine 

w(x.  t)  -  R(h), 

(12) 

with 


R  -  {*,  n.  <?(}■ 

The  Du  hem  inequality  imposes  some  restrictions  on  the  constitutive  relations  (12). 
Assume  that  4  is  a  piecewise  continuously  differentiable  function  with  respect  to  *(x,  t); 


8w 


d<Tu 


(13) 


Then  (4)  becomes 


(,4) 


Since  we  can  choose  arbitrarily  t,  tj  and  &n  for  a  process  and  still  obtain  an  admissible  local 
thermodynamic  process,  it  must  be 


$"°’ 

(15) 

fr“‘n’ 

(16) 

5“°’ 

(17) 

3th 

pT 

p  dw 

(18) 

2.3  Specialization  to  metal  forming  probisms 

By  using  the  concept  of  internal  parameters,  virtually  any  feature  of  material  behavior  can  be  taken  care 
of,  either  macxoscopfcaUy  or  microscopically.  Microscopic  behavior  includes  isotropic  or  anisotropic 
workhardening,  relaxation,  creep,  etc.  and  microscopic  includes,  dislocation  density,  phase  tranformation, 
grain  properties,  etc. 

U  'ortunately.  data  referring  to  this  kind  of  treatment  are  not  easily  available.  Hart  et  of.  [9-11]  have 
dr  pad  an  experimental  program  where  an  exhaustive  number  of  phenomena  are  taken  into  con- 
ration  to  describe  the  nonelastic  deformation  of  metals.  This  type  of  data,  when  fully  available,  will  be 
an  excellent  basis  for  the  present  formulation. 


? 
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For  the  time  being  a  simple  analysis  is  made,  in  *  -ich  the  only  internal  parameters  considered  are  the 
plastic  strain  tensor  and  the  dislocation  density  o.  These  parameters  are  involved  with  a  measurement  of 
workhardening  and  the  change  in  free  energy  that  external  loads  produce  during  deformation  in  a  stable 
non*  recry slalizing  material. 

The  rate  equation  (10)  corresponding  to  the  plastic  strain  tensor  is  simply  borrowed  from  the  mechanical 
analysis  of  rigid  viscoplastic  deformation  [6, 12].  We  recall  that  a  stress  space  yield  surface  expands  with 
the  strain  rate  from  a  static  yield  surface,  and  that  a  Mises-like  law  is  obtained,  provided  we  assume 
isotropy,  isotropic  workhardening,  and  no  further  dependence  on  internal  parameters. 


(19) 


in  (19)  a  is  a  function  of  strain,  strain  rate  and  temperature,  and  is  the  scalar  representation  of  how  the 
yield  surface  behaves.  In  order  to  express  the  dependence  on  strain  rate,  experimental  data  were  found  to 
be  adequately  fitted  by  expressions  of  the  form 

<T“i'[1  +  (y)1'  (20) 


where  Y  is  the  static  yield  stress,  and  Y,  *i  and  y  are  functions  of  strain  and  temperature. 

The  rate  equation  (10)  corresponding  to  the  dislocation  density  is  written  based  on  physical  grounds, 
although  the  values  can  be  accurate  only  to  the  order  of  magnitude.  For  this  reason  we  use  a  simple 
macroscopic  approach,  instead  of  the  elaborate  physical  metallurgy  models.  Fortunately,  as  will  be  seen  in 
the  final  form  of  the  heat  balance  equations,  the  contribution  of  these  terms  is  small.  We  take  then 

a  *  At-  Bat~UkT, 

o(l#)*a»,  (21) 

in  a  simple  attempt  to  say  that  dislocation  density  will  have  both  a  positive  rate  component,  due  to  an 
increase  in  plastic  deformation  (strain  hardening),  and  a  negative  one.  In  terms  of  high  temperature 
deformation,  this  negative  component  can  be  due  to  annihilation  occuning  from  the  diffusion-driven  climb 
of  dislocations.  In  such  a  case,  the  greater  the  number  of  dislocations,  the  higher  the  probability  of 
annihilation. 

The  constitutive  relations  (12)  are  specialized  as  follows. 

A  hot  metalworking  process  does  not  usually  take  very  long,  and  the  temperature  changes  come  from 
plastic  work  and  losses  to  the  environment.  Therefore,  an  isotropic  Fourier  law  for  heat  flux,  with  constant 
thermal  conductivity  k,  is  appropriate 


• 

As  is  done  in  elastic-plastic  theory[13-15],  we  assume  the  following  decomposition 

(22) 

and  take' 

(23) 

(24) 

where  the  specific  heat  (at  constant  volume)  c  is  assumed  constant  in  the  working  temperature  range. 

Now  the  internal  parameters  can  be  used  to  describe  the  energy  which  goes  into  the  structure  of  the 
material.  Since  we  neglect  elastic  deformation,  strain  energy  doesn't  apply,  and  only  the  energy  of 
dislocations  is  involved,  as 


*‘7 

(25) 

From  (23),  (24)  and  (16)  we  have 

*--£(£*")"*?■ 

(23a) 

(16a) 

and  the  dissipation  inequality  (18)  becomes 

alt + Zuk  _  id  i + ife  e-»*r  a  o 
pT  p  p  p 

(26) 

In  (26)  «?,Tj  *0,  therefore. 

Ofjiq  i  vAi  -  vBcte~ukT. 

(26a) 
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This  restriction  is  satisfied  with  experimentally  determined  values  of  v,  A  and  B.  In  fact,  the  r.h.s. 
represents  the  rate  of  energy  spent  in  creating  new  dislocations,  and  it  has  been  shown  experimentally  that 
this  dislocation  energy  ranges  between  S  and  15%  of  the  work  spent  in  plastic  deformation. 

3.  LOADING  AND  BOUNDARY  CONDITIONS 

Usually  temperatures,  temperature  rates  and  stresses  are  not  specified  in  a  process;  instead  body  forces 
and  boundary  tractions,  together  with  heat  supply  both  internal  and  through  the  surface  are  known.  The 
equivalence  is  obtained  through  equilibrium  and  energy  balance  equations. 

As  closed-form  solutions  are  virtually  impossible  to  obtain  in  a  general  problem,  numerical  ap¬ 
proximations  are  sought.  A  finite  element  discretization  is  a  widely  used  and  well-proved  method  to  which 
we  will  refer,  without  entering  into  its  details. 

Assuming  that  no  body  forces  and  no  inertial  effects  are  present,  equilibrium  equations  become,  through 
the  volume 


=■  0. 


(27) 


Writing  equation  (27)  in  a  weak  form 


j  CTiijSu,  dv  =  0. 


Applying  the  divergence  theorem  to  equation  (28),  we  obtain 


<r<i&*n  dv-  J  TM  dSF  =  0, 


(28) 


(29) 


with  imposition  of  boundary  tractions 


T  =  it,,/!,.  (30) 

Boundary  displacements  are  imposed  in  the  form  of  their  rates  (velocities).  By  discretizing,  using  09). 
and  solving  finite  element  equations,  a  velocity  field  is  obtained  and  strain  rates  and  stresses  are 
determined. 

The  energy  balance  equation  (3),  considering  no  external  body  heat  supply,  becomes 

+  k,  TM  -  pet  -  vAi  +  vBae~akT  =0,  (31) 

and,  in  a  weak  form, 

f  k,  TmST  dv  -  f  pctST  dv  +  [  (mfc  -  vAi  +  vBaeakT)ST  do  =  0.  (32) 


By  using  the  divergence  theorem 


£  k,  TiSTi  do  +  £  pctST  do  -  j (oii«ii  -  vA  i  +  vBae'atT)STdv  -  J  q.ST  ds9  «  0, 


(33) 


where 


q.-kj*  (34) 

is  the  boundary  heat  flow. 

Simple  forms  of  boundary  heat  flow  suffice  for  the  present  formulation.  Radiation  heat  flow  is  expressed 
as 


q.'  »<7»(T,‘-T,4).  (34a) 

Convection  heat  flow  is  expressed  as 

q.c  *  h(T,  -  T,).  (34b) 

Temperatures  also  can  be  directly  imposed.  Upon  discretization  and  solution  of  finite  element  equations, 
we  obtain  a  temperature  field  provided  we  know  stresses  and  strain  rates. 

Equations  (29)  and  (33)  are  highly  coupled,  making  a  simultaneous  solution  of  their  finite  element 
counterparts  necessary.  This  process  is  treated  with  time  dimension  analysis. 

4.  TIME  DIMENSION 

In  order  to  describe  a  deformation  process  in  time,  integration  of  rate  equations  (19),  (21)  and  (33)  must 
be  performed  with  respect  to  time.  As  far  as  the  strain  equation  is  concerned,  the  assumption  of  no  inertia) 
effects  requires  that  for  small  intervals  of  time,  velocities  (and  strain  rates)  change  by  very  small  amounts. 
The  strains  at  time  f  +  At  are  obtained  from  their  values  and  rates  at  time  t 


<r»v  “  t,  +  t,Al. 


(35) 


704 


N.  Rebelo  and  S.  Kobayashi 


Similar  treatment  is  made  with  respect  to  updating  displacements,  using  velocities. 

In  standard  finite  element  books,  e.g.  [16],  it  can  be  seen  that  heat  balance  equations  such  as  equation 
(33),  upon  finite  element  discretization  are  reduced  to  the  form 

CT  +  KT  =  Q.  (36) 

The  theory  necessary  to  integrate  these  equations  can  be  found  in  numerical  analysis  books  [17, 18],  It 
suffices  to  say  that  one-step  method  is  used.  Convergence  of  a  scheme  requires  consistency  and  stability. 
Consistency  is  satisfied  by  an  approximation  of  the  type 

=  T,  +  Af •[(  I  -  0)T ,  +  flt,.*] 
where  0  is  a  parameter  varying  between  0  and  I. 

For  unconditional  stability,  0  should  be  greater  than  0-5,  and  a  value  of  0-75  was  chosen.  Selection  of  a 
proper  value  for  0  is  an  important  factor  in  the  situation  where  it  is  desired  that  the  time  step  be  as  large  as 
possible,  provided  that  the  increments  in  strain  are  compatible  with  an  infinitesimal  analysis.  In  the  standard 
test  problem  for  checking  numerical  integration  schemes,  this  value  of  0  gives  good  accuracy  in  a  wide 
range  of  values  of  A t  starting  at  zero. 

An  implementation  of  the  algorithm  above,  considering  T,,u  as  a  primary  dependent  variable[I9],  has 
been  elaborated  here  in  order  to  include  the  coupling  of  the  solution  of  stress  equilibrium  equations  with 
heat  balance,  and  to  include  updating  of  internal  parameters. 

The  sequence  of  calculations  to  be  performed  in  one  time  step  Af  is  as  follows;  (a)  assume  that  the 
initial  temperature  field  and  dislocation  density  distribution  T0  and  at  are  known;  (b)  calculate  initial 
velocity  field  from  equations  (29);  (c)  calculate  initial  temperature  rates  T„,  from  equations  (36);  calculate 
initial  dislocation  density  rates  a0,  from  equations  (21);  (d)  calculate 


(e)  update  displacements  and  strains;  (f)  use  the  old  velocity  field  for  first  approximate  temperature 
calculations 

R  =  Qi'-Ct 

solve 


(*  +  0SfC)*r"“RforT*: 

(g)  calculate  new  velocity  held  from  equations  (29):  (h)  calculate  second  values  for  temperatures 

R-QS'-Ci 

solve 

(k+^c)*T1u  “RforT2' 


use  a  similar  procedure  to  calculate  dislocation  density;  (i)  iterate  until  convergence:  (j)  calculate  new 
temperature  rate 


tw 


and  new  dislocation  density  rate. 

It  should  be  noted  that  integration  of  equations  (21)  follows  the  same  procedure  described  above,  since 
equations  (21)  and  (36)  are  of  the  same  type. 


5.  CONCLUDING  REMARKS 

In  order  to  simplify  the  description  of  the  algorithm  for  the  coupled  analysis,  the 
solution  for  the  velocity  fields  has  not  been  detailed.  The  finite  element  equations 
obtained  from  (29)  and  ( 19)  are  nonlinear,  and  an  iterative  scheme  based  on  the  Raphson 
method  is  used.  Such  a  scheme,  as  well  as  convergence  criteria,  has  been  described  at 
length  in  previous  papers  [2, 3, 5], 

Dislocation  density  values  used  are  averages  for  each  element.  This  procedure  is 
justified  by  the  large  uncertainty  in  data,  and  by  the  moderate  influence  of  these 
values  in  the  rest  of  the  calculations. 

We  have  presented  the  theoretical  background  necessary  for  the  solution  of 
metalworking  problems  in  which  temperature  variations  as  well  as  a  heat  transfer 
analysis  are  important.  The  concept  of  internal  parameters,  and  the  way  they  are 
imbedded  in  the  solution  algorithm,  opens  good  perspectives  for  better  descriptions  of 
material  behavior,  and  of  physical  phenomena  concurrent  with  plastic  deformation. 
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Sammary — The  finite  element  formulations  for  coupled  analysis  of  deformation  and  heat 
transfer  are  applied  to  the  problems  of  solid  cylinder  compression  and  ring  compression.  The 
temperature  distributions  in  solid  cylinder  compression  are  compared  with  experimental  data, 
and  the  computed  results  of  ring  compression  of  aluminum  at  elevated  temperature  are 
compared  with  experimental  observations  as  well  as  with  calculated  results  not  including  the 
temperature  effects.  It  is  demonstrated  that  the  coupled  analysis  accurately  simulates  these 
processes. 


NOTATION 

Y  static  yield  stress 
«  strain  rate 
y,  n  material  parameters 
Yn  static  yield  stress  at  reference  temperature  T 
Q  material  parameter 
k  Boltzman  constant 
T  absolute  temperature 
f  heat  capacity  matrix 
£  heat  conduction  matrix 
Q  heat  flux  vector 
T  vector  of  nodal  point  temperatures 
T  vector  of  nodal  point  temperature  rates 
N  interpolation  function 
ai  rate  of  plastic  work 

vAi  -  vBa  e'“*T  rate  of  change  of  dislocation  energy 
A,  B  material  constants 
a  dislocation  density 
D  energy  for  self  diffusion 
v  energy  per  unit  length  of  dislocation 
T.  environment  temperature 
T,  surface  temperature 
T*  die  temperature 
T.  workpiece  temperature 
<r  Boltzman  constant 
<  radiation  coefficient 
h  heat  convection  coefficient 
fiu  heat  transfer  coefficient  from  lubricant 
4/  surface  heat  generation  rate  due  to  friction 
/  friction  force 

c,  relative  velocity  between  the  die  and  the  workpiece 

H  workpiece  height  in  mm 

*.  heat  conductivity— workpiece 

K  heat  conductivity— die 

c»  heat  capacity— workpiece 

c4  heat  capacity— die 

p  mass  density 

n  number  of  strokes  per  minute  (press) 

S  stroke  of  the  press 
W  distance  from  lower  position 


1.  INTRODUCTION 

The  theory  necessary  for  developing  computer  simulations  of  transient  viscoplastic 
metal  deformation  coupled  with  a  heat  transfer  analysis  has  been  presented  in  the 
companion  paper  [1], 
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When  deformation  takes  place  at  high  temperatures,  material  properties  caa  vary 
considerably  with  temperature.  During  a  metalworking  process,  large  aonhomogeweities 
in  deformation  and  consequently  in  heat  generation,  usually  occur.  Moreover,  especially 
if  the  dies  are  at  a  considerably  lower  temperature  than  the  workpiece,  the  heat  losses  by 
conduction  to  the  dies  and  by  radiation  and  convection  to  the  environment  contribute  to 
the  existence  of  severe  temperature  gradients.  Thus  including  temperature  effects  in  the 
analysis  of  metal  forming  problems  is  very  important. 

Furthermore,  at  elevated  temperatures  plastic  deformation  can  induce  phase 
transformations  and  alterations  in  grain  structures  which  in  turn  will  modify  the  flow 
resistance  of  the  material  as  well  as  other  mechanical  properties.  Including  these 
metallurgical  effects  in  the  process  analysis  contributes  significantly  to  the  under¬ 
standing  of  material  behavior  under  plastic  deformation. 

The  importance  of  temperature  calculations  during  a  metal  forming  process  has 
been  recognized  for  a  long  time.  Until  recently  the  majority  of  the  work  has  been 
based  on  procedures  that  uncouple  the  problem  of  heat  transfer  from  the  metal 
deformation  problem. 

Several .  esearch  workers  have  used  the  approach  of  determining  the  flow  charac¬ 
teristics  of  the  problem  either  experimentally  or  by  calculations,  and  then  of  using 
these  calculation  for  the  calculation  of  heat  generation.  Among  these  is  the  work  of 
Johnson  and  Kudo [2]  on  extrusion,  and  Tay,  Stevenson  and  Dans(3]  on  machining. 
Another  approach  uses  Bishop’s  [4]  numerical  method,  in  which  heat  generation  and 
transportation  is  considered  to  occur  instantaneously  for  each  time  step,  with  con¬ 
duction  taking  place  during  the  time  step.  The  new  temperature  distribution  is  used  to 
calculate  flow  characteristics,  which  in  turn  will  originate  the  instantaneous  heat 
generation  and  transportation  for  the  next  step.  This  approach  is  represented  by  the 
work  of  Aitan  and  Kobayashi[$]  on  extrusion,  Lahoti  and  Altan[6]  on  compression 
and  torsion,  and  Nagpal,  Lahoti  and  Aitan[7]  on  forging.  Usually  the  temperature 
calculations  are  done  by  finite  differences,  or  by  finite  elements,  and  the  upper  bound 
technique  is  the  most  common  for  determining  flow  patterns,  if  not  experimental. 

As  far  as  we  know,  a  paper  by  Oden,  Bhandary,  Yagawa  and  Chung[8]  presents  the 
first  coupled  analysis  of  deformation  and  heat  transfer,  although  the  application  is  to  a 
3-dimensional  rectangular  bar,  which  was  constrained  on  one  direction  and  heated  at  a 
corner.  The  formulation  is  for  elastoviscoplastic  materials,  and  uses  the  same  rational 
thermodynamics  background  we  have  used  in  the  present  work. 

Zienkiewicz,  Onate  and  Heinrich  [9]  have  made  a  coupled  thermal  analysis  in 
extrusion.  Because  it  is  a  steady  state  situation,  the  heat  balance  equations  and 
equilibrium  equations,  using  rate  sensitive  materials,  are  solved  simultaneously. 
Previously  Zienkiewicz,  Jain  and  Onate  [10]  had  proposed  to  solve  the  same  problem 
by  separating  the  equations,  and  iterating  the  results  between  the  two  sets  of  solutions. 

The  purpose  of  this  paper  is  to  show  the  results  of  calculations  by  applying  the 
formulation  developed  in  [1]  to  metal  forming  problems,  which  demonstrate  tem¬ 
perature  effects  clearly,  and  to  show  the  potential  for  further  refinements.  A  brief 
discussion  of  the  basic  finite  element  equations  is  presented,  in  which  some  details  of 
the  treatment  of  boundary  conditions  crucial  for  the  algorithm  are  given. 

The  compression  of  a  solid  cylinder  of  AISI  1013  is  analyzed  and  compared  with 
experimental  data  obtained  by  Pohi[lI].  Also,  the  computed  results  of  the  forging  of 
an  aluminum  1 100  ring  at  400°C  are  compared  with  the  one  studied  in  [7]  as  well  as  with 
calculations  hot  including  the  temperature  effects  [12]. 

2.  BASIC  EQUATIONS 

The  finite  element  equations  resulting  from  discretizing  weak  forms  of  equilibrium  equations  and  heat 
balance  equations  are  described  in  the  literature.  We  discuss  here  only  the  details  that  are  important  in  the 
implementation  of  the  algorithm  presented  in  the  companion  paper. 

As  far  as  equilibrium  equations  are  concerned,  the  only  modification  of  the  analysis  of  Oh.  Rebeto  and 
Kobayashi(l2]  is  the  temperature  dependence  of  the  flow  stress.  Recall  that  experimental  data  on  flow 
stress  has  been  fitted  by  the  expression  of  the  form 


lit  of  vitcopiartic  deformation — II 


709 


and  that  Y.  a,  y  are  feartinnr  of  strain  tod  temperature.  For  the  range  of  temperature  .nations  in  the 
metal  forming  problen,  and  within  the  accuracy  of  dau  available,  equation  (1)  can  be  e  ‘ssed  as 

* ['♦(*)  ]•  <2> 

Y*.  a  and  y.  functions  of  strain 

The  friction  factor  is  dadoed  at  the  ratio  batwooa  the  fncttoaal  shear  stress  and  the  rate  and 
ToHmm  bnlancraqiiationa,  see  recall  again  that  there  are  of  the  form  (see  E.G.[13)); 

ct-rr-Q.  (3) 

where  the  heat  flux  vector  Q  has  several  components  (tee  (11); 

Q-Jfw#- vM  +  vflo  r‘,*lT')N dV  +  J  T/)N dS, 

♦  J  *(T.  -  T,)N  dS,  +  £  hwfli  ~  r.)N  dS,  «> 

+  f  pfNdSi. 

The  first  component  is  the  contribution  of  the  net  heat  federated  inside  the  deforming  body,  the  second 
one  the  contribution  of  the  heat  radiated  into  the  body,  and  the  third,  the  contribution  of  the  heat 
transferred  by  convection.  The  fourth  term  represents  the  contribution  of  the  heat  transferred  between 
workpiece  and  die  through  their  interface  (with  or  without  lubricant),  and  the  last  term  is  the  contribution  of 
the  heat  peneraled  by  friction  between  die  and  workpiece. 

The  following  remarks  pertain  to  the  implementation  of  the  algorithm  previously  described. 

We  treat  the  workpiece  and  the  die  separately.  In  this  way,  we  greatly  reduce  the  number  of  equations 
to  be  solved  simultaneously,  which  in  turn  reduces  the  cost  of  the  solution.  In  the  die  there  is  no  interna] 
heat  feneration  and  no  deformation  calculations  are  necessary  (rigid  dies).  Because  of  the  heat  transfer  at 
the  interface,  the  die  and  the  workpiece  equations  contain  temperatures  of  each  other,  and  the  solution  is 
found  iteratively.  For  two  reasons,  this  brings  no  increase  in  cost:  first,  iterations  with  respect  to  the 
velocity  have  to  be  done  anyway;  and  second,  matrices  do  not  have  to  be  recalculated  at  each  iteration, 
as  only  the  heat-flux  vector  changes. 

In  a  numerical  testing  of  heat  transfer  calculations,  it  was  found  that  lubricant  heat  transfer  coefficients 
at  least  one  order  of  magnitude  higher  than  the  real  ones  were  necessary  to  produce  instability. 

The  heat  generated  through  friction,  and  calculated  by 

9/  -  I/I  •  kl. 

was  evenly  distributed  between  the  die  and  the  workpiece. 

The  die  had  to  be  remeshed  at  each  time  step,  in  order  to  make  the  nodes  at  the  interface  match.  Again, 
numerical  testings  were  made  with  similar  mesh  changes  in  a  heat  transfer  problem  whose  exact  solution  is 
known,  and  the  error  in  temperatures  so  introduced  was  very  small. 

The  temperature  values  used  in  the  internal  heat  generation  and  radiation  terms  of  heat  flux  are  at  all 
times  those  values  from  the  previous  iterations. 

Dislocation  density  calculations,  made  only  in  the  workpiece,  represent  a  minute  fraction  of  operations, 
because  no  coupling  was  considered  from  element  to  element.  For  each  element  the  scalar  version  of  the 
time  integration  scheme  described  for  temperatures  was  used. 

3.  COMPRESSION  OF  A  CYLINDER  OF  AJSI  1013  AT  ROOM  TEMPERATURE 

Pohl(U]  conducted  temperature  measurements  in  order  to  test  his  uncoupled  analysis,  in  which 
approximate  stream  functions  were  used  for  the  deformation,  and  finite  differences  for  the  heat  balance.  A 
solid  cylinder  of  a  carbon  steel  AISI  1013  was  compressed  between  flat  dies  at  room  temperature. 
Thermocouples  were  inserted  in  the  cylinder  at  different  locations.  Upon  deformation,  their  readings 
indicated  the  temperature  increases  due  to  the  heat  generated,  but  no  assessment  was  made  with  respect  to 
the  accuracy  obtained  in  such  measurements.  Nevertheless,  this  is  one  of  few  cases  in  the  literature  where 
actual  temperature  distributions  were  reported. 

Fig.  1  shows  the  dimensions  and  locations  of  measuring  points. 

Although  the  values  of  some  parameters  must  be  “guessed”  in  computation,  the  results  of  our  analysis 
are  compared  with  those  measurements. 

The  conditions  used  in  the  computation  are  as  follows. 

The  deformation  took  place  at  initial  room  temperanne,  up  to  a  height  reduction  of  1/3.  The  finite  element 
grid  was  composed  of  132  four  node  quadilateral  elements  in  the  workpiece,  and  1 19  in  the  die.  Because  of 
symmetry,  only  one  fourth  of  the  cylinder  needed  to  be  analyzed.  The  friction  factor  was  taken  as  0-63. 

In  order  to  simulate  a  mechanical  press,  the  die  velocity  was  changed  at  each  time  step  according  to  the 
expression 

v  m  12  *  V(H-20)  mm/s. 

Each  step  corresponded  to  Iffi  in  height  reduction,  which  was  equivalent  to  time  steps  of  up  to  0*03  sec. 

In  this  case  the  flow  stress  was  taken  independently  of  strain  rate  and  temperature.  This  is  justifiable  for 
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Fig.  I.  Specimen  dimensions.  Die  dimensions.  Locations  of  temperature  measurements/ll]. 

steel  being  deformed  at  room  temperature.  The  flow  stress  values,  as  given  by  Pohi,  were 

<r»  275  N/mm  e  <002312 

<r-  722*  r**2  N/mm  f  a  0  025 12. 

The  heat  transfer  characteristics,  other  than  the  thermal  conductivity  and  heat  capacity  of  AISI  1015 
measured  by  Pohl,  are  taken  from  handbooks; 

k.  «*  36N/s°K 
pc.  *  3-77  N/mm2  “K 
t,«19N/s*K 
pct  *  3-77  N/mm1  “K 
era- 85*10-“ N/mms  fK/ 
h«  0  00295  N/s  mm  “K 
h)l*«4N/Smm,K. 

The  dimensions  of  the  die  are  such  that  at  the  outside  boundaries  a  constant  temperature  is  imposed. 

The  material  was  given  an  initial  dislocation  density  of  lOVmm2.  and  constants  A,  B,  v  and  D  were 
chosen  in  such  a  way  that  practically  no  annihilation  occurs  (room  temperature),  and  the  density  grows 
with  strain  to  an  order  of  lOtymm2  at  a  strain  of  about  50% 

A  -  125*  lOVmm2 
B«I55*IO»/s 
Dfk-  17000  *K 
v  «  317*  10^"  N  mm/mm. 

Fig.  2  shows  the  initial  finite  element  layout,  and  Fig.  3  shows  how  it  deformed  after  1/3  reduction  in 
height.  The  barreling  and  the  nearly  sticking  situation  along  the  interface  should  be  noted. 

In  Fig.  4  the  effective  strain  distribution  is  shown,  which  points  out  the  nonuniformity  in  deformation. 
The  most  deformed  zone  is  the  outside  corner,  and  shear  band  formation  is  clearly  seen.  As  expected,  the 
zone  under  the  die  shows  little  deformation. 

Temperature  contours  in  both  workpiece  and  die  are  presented  in  Fig.  5.  In  the  workpiece  the 
temperature  distribution  has  a  similar  pattern  to  that  of  effective  strain  distribution.  Since  the  initial 
temperatures  were  the  same  for  die  and  workpiece,  and  equal  to  room  temperature,  internal  heat  generation 
is  clearly  the  prevailing  factor  in  modifying  temperatures.  The  temperature  distribution  in  the  die  shows 
well  that  the  most  deformed  zone  was  the  largest  source  of  heat. 

The  temperature  values  measured  by  Pohl  are  compared  with  the  calculated  values  in  Figs.  6-8. 

The  agreements  are  excellent  for  the  inner  points,  as  seen  in  Figs.  6  and  7.  For  the  most  outer  point!, 
however,  the  computed  results  indicate  that  the  temperature  difference  for  the  three  points  is  minor,  while 
experiments  show  more  difference.  Nevertheless,  this  difference  does  not  seem  to  negate  the  analysis 
technique,  as  the  heat  conduction  part  of  the  program  has  been  tested  with  accurate  results  and  the 
accuracy  of  the  finite  element  deformation  analysis  has  been  well  established. 

The  discrepancies  seen  in  Fig.  8  can  be  attributed  to  possible  inappropriate  material  constants  and 
experimental  inaccuracies. 
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Flo.  4.  Effective  strain  distribution  after  1/3  reduction  in  height 
FKJ.  3.  Temperature  distribution  after  1/3  reduction  in  height 
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Flo.  6.  Comparison  of  lempira  two  distribution  between  theory  (— )  and  experiments  (x.O,  +). 
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Fio.  7.  Comparison  of  temperature  distribution  between  theory  (-)  and  experiment  (*.+). 
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Fig.  8.  Comparison  of  temperature  distribution  between  theory  ( — )  and  experiment  (x,  O,  +). 

4.  COMPRESSION  OF  A  RING  OF  ALUMINUM  1100  AT  HIGH  TEMPERATURE 

Complete  finite  element  solutions  for  ring  compression  of  rigid  plastic  materials  have  been  obtained  by 
Chen  and  Kobayashi[l4].  In  [12]  they  were  extended  to  rate-sensitive  materials.  The  purpose  of  analyzing 
the  problem  with  temperature  effects  is  twofold;  to  compare  calculations  with  some  experimental 
measurements  of  Nagpal,  Lahoti  and  Altan[7]  and  to  compare  them  with  our  previous  results  without 
temperature  effects.  The  success  of  the  analysis  makes  it  possible  to  extend  the  ring  compression  test  to 
quantitative  lubricant  evaluation  in  the  hot  forming  range. 

The  rings  tested  have  a  o.d. ; i.d. ; height  ratio  of  6:3:2.  and  an  outer  diameter  of  54  mm. 

As  both  dies  and  workpiece  had  their  surfaces  cleaned,  and  no  lubricant  was  used,  the  friction  factor 
was  taken  as  10.  The  dies  were  initially  at  room  temperature,  and  the  rings  were  heated  to  427*C  (800° F). 

Total  reduction  in  height  was  50%,  and  the  calculations  were  performed  in  steps  of  1%  deformation 
each.  T.x:  equivalent  time  steps  did  not  exceed  0  003  sec.  Assuming  symmetry,  one  quarter  of  the  ring  was 
analyzed  and  104  4-node  quadrilateral  elements  were  placed  in  the  ring  and  150  in  the  die.  Fig.  9  shows 
the  initial  grid  layout  as  well  as  the  dimensions. 

In  order  to  simulate  a  mechanical  press,  the  instantaneous  velocities  of  the  die  were  obtained  from  the 
expression  [15]; 


v  -  0- 1  •  n  •  V(S  •  W  -  W*)  mm/s 
n  *90 
5 -254  mm 
W  *9  mm  initially. 

Flow  stress  expressions  for  Aluminum  1100  were  derived  from  experimental  data  of  Lindholm[16]  and 
Altan  and  Boulger[17],  The  expressions  in  Table  I  correspond  to  a  reference  temperature  of  400*C  to  which 
the  following  factors  in  equation  (2)  were  applied 

ew*T*  -  3-55 
Qlk  «  852-6855*K. 


Table  1.  Stress- strain  strain -rate  properties  of  A1  1100  at 
400*C  USED  IN  CALCULATIONS 
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Flo.  9.  Initial  mesh  for  Al  1100  ring  compression. 
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The  heat  transfer  characteristics  were  taken  from  standard  handbooks. 

fc.  •  242  N/s  *K 
pc. -2-4255  N/mm2*K 
kt  -  19  N/s  *K 
pc<- 3-77  N/mm3*K 
<re-85*IO-,>N/mms*K< 

A -0  007  N/s  mm  *K 
Am -35  N/s  mm  gK. 

As  Fig.  9  shows,  the  outer  boundaries  of  the  die  are  sufficiently  far  from  the  interface  so  that  a  constant 
temperature  is  imposed  on  them. 

The  dislocation  data  utilized  were  prepared  under  two  basic  assumptions.  First,  if  there  were  no 
annihilation,  dislocation  density  would  grow  from,  say,  lOVmm2  in  the  annealed  state  to  HP/mm2  at  an 
effective  strain  of  50%.  Second,  because  of  annihilation  at  high  temperatures  and  in  a  relatively  fast 
process,  dislocation  density  will  not  grow  above  IO*/mm2.  The  values  used  are 

A»2»I<P/mm3 
ff  -  I  •  I  •  10,4/s 
m- 17000  *K 
v  -  105  *  10"  N  mm/mm. 

in  deformation  calculations  an  average  of  12  iterations  were  needed  for  convergence,  and  up  to  3  more 
for  temperatures.  Whenever  a  node  on  the  surface  touched  the  die,  convergence  was  much  much  slower 
and  up  to  30  iterations  were  then  needed. 

Fig.  10  compares  the  decrease  in  inner  diameter  vs  reduction  in  height  obtained  in  the  present 
calculations  with  the  experimental  point  of  Nagpal  el  aL  and  with  the  previous  solutions  not  including 
temperature  effects.  When  the  effect  of  temperature  is  included,  the  effect  of  chilling  by  the  dies  causes 
relatively  small  sliding  along  the  interface  and  results  in  more  decrease  in  inner  diameter  than  the  case 
without  temperature  effect 

Fig.  II  shows  the  ring  profiles  after  48-5%  reduction.  The  agreement  is  excellent  between  computed 
and  experimental  profiles. 

Grid  distortions,  indicating  the  How  of  material  after  20  and  50%  reduction,  are  presented  in  Figs.  12(a) 
and  13(a).  The  neutral  zone  is  distinctively  traceable,  and  the  large  amounts  of  folding  are  evident.  Then- 
equivalent  effective  strain  distributions,  clearly  showing  the  nonuniform  deformation,  are  given  in  Figs. 
14(a)  and  15(a).  The  most  deformed  areas  are  the  original  corners  in  contact  with  the  die,  and  the  least 
deforming  regions  are  noticeable  near  the  neutral  axis.  A  comparison  with  equivalent  figures  (12(b)—  15(b)) 
for  the  temperature-independent  calculations  in  similar  conditions  indicates  a  more  severe  distortion  in  the 
temperature  dependent  case,  and  demonstrates  how  the  die-c hilling  effect  concentrates  the  deformation  on 
the  equatorial  regions,  which  are  the  hottest  and  least  resistant  to  flow. 

A  characteristic  behavior  of  calenlated  dislocation  densities  is  that  they  increase  sharply  during  the  first 
5%  of  deformation,  to  value*  from  MP/mm2  up  to  5*  lO’/mm2  At  this  point  the  elements  at  higher 
temperatures  begin  to  have  negative  rates,  while  the  ones  near  the  die  keep  increasing  their  densities.  At 
25%  deformation  the  hipest  value  of  1-2*  HP/mm2  is  reached,  in  the  upper  right  corner  element.  By  that 
time  the  great  majority  of  elements  have  densities  between  I  and  3  •  lO’/mm3.  although  some  are  increasing 
and  others  decreasing  When  46%  deformation  is  reached.  aR  elements  start  to  have  negative  rates  of 
dislocation  density.  This  means  that  the  density  values  reached  are  high  enough  to  induce  strong  annihilation, 
even  in  the  coUsr  zones  near  the  die.  The  average  values  at  this  stage  are  in  the  HP/mm2  order  of  magnitude  for 
the  hottar  zones,  and  the  lO’/mm2  order  of  magnitude  for  the  colder  zones. 


Flo.  10.  Chengs  or  ring  inner  diameter  as  functions  of  reduction  in  height. 
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Fig.  II.  Deformed  ring  profiles. 
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Fig.  12.  Grid  distortions  of  Al  1100  ring  at  20%  reduction  in  height  (a)  with  temperature  effect 
and  (b)  without  temperature  effect. 


Flo.  13.  Grid  distortions  of  Al  1100  ring  at  50%  reduction  in  height  (a)  with  temperature 
effect  and  (b)  without  temperature  effect. 


FW.  14.  Strain  distributions  of  A I  1100  ring  at  20%  reduction  in  height  (a|  with  temperature 
effect  and  (b)  without  temperature  effect. 
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Fig.  15.  Strain  distribution  of  At  1100  ring  at  50%  reduction  in  height  (a)  with  temperature  effect 
and  (b)  without  temperature  effect. 

We  chose  the  material  constants  for  the  dislocation  density  calculations,  having  in  mind  that  in  overall 
deformation  the  energy  spent  with  the  creation  of  dislocations  is  around  10%  of  the  total  deformation 
energy.  The  results  discussed  above  show  that  the  former  energy  is  mostly  spent  in  the  early  stages  of 
deformation,  and  that  towards  the  end,  annihilation  is  helping  plastic  work  to  supply  heat  generation 
throughout  the  body. 

In  Figs.  16  and  17  temperature  contours  in  both  ring  and  die  are  presented.  The  temperature  field  is 
dominated  by  the  heat  flow  from  the  workpiece  to  the  die.  and  only  near  the  equator  temperatures  rise 
slightly.  Temperature  differences  across  the  workpiece  die  interface  decrease  from  407  K  initially,  to 
around  220  K  at  the  end  of  deformation. 

An  appraisal  of  influence  of  parameters  in  the  temperature  calculations  was  done  by  roughly  calculating 
heat  fluxes  in  I  radian  section  for  a  typical  step  of  deformation,  and  comparing  them  to  heat  generation; 

radiated  heat  -  350  Nnun/i 
converted  heat  -  500  Nmm/s 
beat  lost  to  die  ~  6500000  Nmm/s 
heat  generated  —  1000000  Nmm/s 

this  indicates  that  heat  feneration  and  conduction  to  the  die  govern  the  whole  heat  balance  process. 

In  Fig.  18  the  load  displacement  curve  is  compared  to  the  previously  calculated  one  not  including 
temperature  effects.  During  the  initial  stages  of  deformation  the  values  are  similar,  due  to  the  fart  that  in 


Fig.  16.  Temperature  contours  at  20%  reduction  in  height  in  Al  1 100  ring  compression. 


Fto.  17.  Temperature  contours  at  50%  reduction  in  height  in  Al  1 100  ring  compression. 
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Fig.  18.  Load  displacement  curves  for  At  1100  ring  compression. 


the  temperature  independent  calculations  the  flow  stress  used  referred  to  a  temperature  similar  to  the  initial 
one  in  the  temperature  dependent  case.  As  deformation  proceeds,  the  zones  near  the  die  cool  down  and  the 
equatorial  zones  heat  up.  Since  deformation  occurs  mainly  in  the  latter,  the  lower  flow  stress  here  induces  a 
lower  total  load  necessary  for  deformation. 

The  decrease  in  load  towards  the  end  is  a  result  of  the  rapid  decrease  in  velocity  of  the  die,  and  the 
subsequent  strain  rate  effects. 


5.  CONCLUDING  REMARKS 

The  capabilities  of  a  coupled  thermoviscoplastic  analysis  of  non-steady-state  metal 
forming  processes  have  been  demonstrated  in  this  paper.  The  finite  element  analysis 
of  deformation  and  heat  transfer  is  a  significant  extension  to  the  analysis  of  warm  and 
hot  forming  processes.  It  produces  accurate  simulations  of  the  processes  and  provides 
insight  into  the  governing  parameters  in  the  process. 

A  higher 'level  of  sophistication  in' the  analysis  is  invariably  associated  with  a 
higher  level  of  sophistication  in  the  input  data  necessary.  At  this  stage  much  work  is 
still  needed  in  order  to  utilize  all  the  capabilities  of  the  method.  It  is  hoped  that  in  the 
future  the  important  internal  parameters  will  be  screened  from  the  superfluous,  and 
programs  to  quantify  their  importance  will  be  set. 
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#3. 

An  Investigation  of  the  Shell 
Nosing  Process  by  the  Finite- 
Element  Method.  Part  1:  Nosing 
at  Room  Temperature  (Cold 


The  metal-forming  process  of  shell  nosing  at  room  temperature  was  analyzed  by  the 
finite-element  method.  The  strain-rate  effects  on  materials  properties  were  included 
in  the  analysis.  In  cold  nosing  simulations;  the  nine-node  quadrilateral  elements 
with  quadratic  velocity  distribution  were  used  for  the  workpiece.  The  treatment  of  a 
moving  boundary  in  the  analysis  of  nosing  is  discussed  and  successfully  im¬ 
plemented  in  the  finite-element  program.  FEM  simulations  of  105 -mm  dia.  shells  of 
AISI 1018  steel  and  aluminum  2024  were  performed  and  solutions  were  obtained  in 
terms  of  load-displacement  curves,  thickness  distribution,  elongation,  and  strain 
distributions.  Comparisons  with  experimental  data  show  very  good  agreement. 


Introduction 

Shell  nosing  refers  to  the  process  of  forming  an  ogive  nose 
at  the  end  of  a  tubular  part  by  pressing  the  tube  into  a  suitably 
formed  die  (Fig.  1).  The  flow  of  metal  in  nosing  is  very 
complex  and  only  a  slight  variation  in  friction  or  temperature 
conditions  may  result  in  a  misformed  shell  due  to  improper 
metal  flow.  Understanding  the  deformation  mechanics  of  a 
metal-forming  process  is  an  essential  and  important  step 
towards  the  control  of  metal  flow,  and,  subsequently,  the 
design  of  the  forming  process.  Experimental  characterizations 
of  process  parameters  involve  trial  and  error,  and  it  is  dif¬ 
ficult  to  separate  the  influence  of  individual  parameters  upon 
each  other.  It  is  desirable  and  necessary  to  use  an  inexpensive 
analytical  simulation  method  in  order  to  improve  industrial 
productivity. 

Successful  efforts  have  been  carried  out  by  Kobayashi  and 
co-workers  (1-4]  using  the  finite-element  method  in  analyzing 
various  metal-working  processes,  such  as  compression, 
heading,  piercing,  and  extrusion.  The  objective  of  the  present 
investigation  is  to  develop  the  finite-element  simulation 
technique,  including  flow  stress  dependence  on  strain  rates, 
for  the  shell  nosing  process. 

The  flow  of  metal  in  the  nosing  process  is  affected  by  the 
characteristics  of  the  lubricant  and  the  speed  of  the  nosing 
operation.  If  the  lubrication  is  not  sufficient  and  if  the  nosing 
process  is  slow,  the  lubricant  is  squeezed  out  and  the  friction 
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between  shell  and  die  tends  to  shorten  the  shell.  On  the  other 
hand,  if  the  lubrication  is  good  and  if  the  nosing  proceeds 
rapidly,  the  lubricant  is  squeezed  out  to  a  lesser  extent. 
Therefore,  the  friction  is  reduced  and  the  shell  is  lengthened. 

Small  caliber  shells  (up  to  105  mm)  can  be  cold  nosed,  while 
larger  shells  (155  mm  and  up)  are  hot  nosed.  In  cold  nosing, 
the  entire  shell  is  initially  maintained  at  a  constant  tem¬ 
perature  and  the  flow  stress  is  uniform.  The  deformation 
work-hardens  the  shell  material  nonuniformly,  and  the 
resistance  to  further  deformation  by  various  parts  of  the  shell 
is  unequal.  The  elastic  spring  back  is  very  small  compared  to 
plastic  strain;  therefore,  the  dies  are  given  the  shape  of  the 
rough-finished  nose.  The  slight  elastic  back  spring  is  removed 
in  the  process  of  finish  turning. 

In  1943  Nadai  [5]  conducted,  for  the  ASME  Special 
Research  Committee  on  Forging  of  Steel  Shells,  an  extensive 
theoretical  investigation  on  the  forces  required  in  nosing  of 
shells.  This  work  is  basically  an  extension  of  the  theory  of 
curved  shells  to  cover  cases  in  which  the  metal  is  in  a  plastic 
state  of  equilibrium.  The  approximation  is  based  upon 
constant  coefficient  of  friction  and  ignores  variations  in  the 
flow  stress  of  the  shell  material.  Onat  and  Prager  [6]  extended 
Nadai’s  work  to  include  the  changes  in  the  shell  wall  thickness 
due  to  the  nosing  operation.  They  investigated  the  influence 
of  these  changes  on  the  stresses  in  nosing.  Cruden  and 
Thompson  [7]  conducted  an  experimental  study  of  the  nosing 
process  to  establish  the  limitations  of  the  process  and  to  assess 
-the  effects  of  the  various  process  parameters.  Recently,  the 
metal  flow  in  nosing  of  shells  was  analyzed,  using  Hill’s 
general  method  of  analysis,  by  Lahoti  and  Altan  (8]. 

In  the  present  study  of  cold  nosing  operation,  the  rigid- 
viscoplastic  finite-element  method  is  used  for  the  analysis, 
with  reference  to  the  effects  of  friction.  Simulation  of  105- 
mm  shell  nosing  for  A1SI  1018  steel  is  performed,  and 
solutions  are  compared  with  available  experimental  data  from 
Battelle  Columbus  Laboratories  [9],  The  experiments  of  shell 
nosing  of  aluminum  alloy  2024  were  conducted  and  an  ex¬ 
tensive  comparison  was  made  between  the  computer 
simulation  and  experiments. 

Finite-Element  Method 


of  the  nodal  point  velocity  vector  v  and  the  element  Lagrange 
multiplier  X,  the  system  of  equations  is  obtained  in  the  form 
of 

f(v,X)-0  (3) 

Equations  (3)  are  nonlinear  in  v  and  linear  in  X.  To  solve  these 
nonlinear  equations,  the  Newton-Raphson  iterative  method 
was  used.  Expanding  equations  (3)  into  Taylor’s  series  and 
maintaining  only  the  first  two  terms, 

f(v,"X)+  (v„,X)Av„  -0  (4) 

which  results  in  the  system  of  equation  linear  in  Av  and  X.  In 
equation  (4),  v„  represents  the  input  values  from  the  previous 
solution  and  Av„  represents  the  correction  values  obtained 
from  the  present  iteration.  After  each  iteration,  the  values  of 
v  are  updated  as  v„+!  *  v„  +  aAv„ ,  where  a  is  a  deceleration 
coefficient.  For  solution  convergence,  two  measures,  namely, 
the  error  norm  of  the  velocities  and  the  norm  of  the  residual 
equation  (3),  are  used  as  criteria. 

The  constitutive  equation  (2)  becomes  undetermined  at  zero 
strain  rate.  Whenever  some  regions  in  a  nonuniform 
deformation  process  become  rigid,  the  system  of  equations 
becomes  very  ill  conditioned.  To  avoid  this,  a  strain  rate 
offset  £q  is  introduced  such  that  the  stress  drops  linearly  to 
zero  when  the  strain  state  is  below  This  .offset  value  is 
chosen  only  for  numerical  reasons.  Values  of  two  orders  of 
magnitude  smaller  than  an  average  strain  rate  in  the  defor¬ 
mation  have  given  successful  results  (11]. 


General  theories  about  the  finite-element  method  have  been 
well  described  in  many  books.  Therefore,  only  a  brief 
description  of  the  method  that  is  necessary  for  the  discussion 
of  the  features  involved  is  given  here.  Detail  formulations  and 
computer  programs  for  the  analysis  can  be  found  in  reference 
[10]. 

For  an  incompressible,  rigid-viscoplastic  material,  the 
variational  formulation  in  boundary  value  problems  is  given 
by 

j  ^  aljUijdV—  js  F,5ufdSF  +  j  ^  X5 itldV 


Flg.2  Math  system  Mid  boundary  condition 


in 


In  equation  (1),  o,y  is  the  deviator  stress,  e<,  is  the  strain ‘rate, 
F,  is  the  traction  prescribed  over  the  surface  Sf,  u,  is  the 
velocity,  and  X  is  the  Lagrange  multiplier. 

Constitutive  equations  for  viscoplastic  materials  are  ex¬ 
pressed  by 


,  2  » 
a‘j  ■  y  j 


(2) 


where  »  is  the  effective  stress  which  is  a  function  of  strain  and 
strain  rate,  and  i  is  the  effective  strain  rate. 

Discretization  of  equation  (1)  is  performed  by  using  nine- 
node  isoparametric  elements  for  the  axisymmetric  defor¬ 
mation  problem.  Applying  the  arbitrariness  of  the  variations 
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Tabic  1  Dimension  of  lOS-mm  ibdl 


Outside  diameter  at  the  nose  base  108  mm 

Radius  of  curvature  of  the  nose  ogive  650  mm 

Initial  wall  thickness  of  the  preform  12  mm 

Axial  length  of  the  nose  portion  ISO  mm 

Distance  of  the  ogive  center  from  nose  base  25  mm 


Computational  Procedures 

In  applying  the  finite-element  method  to  shell  nosing,  some 
thought  should  be  given  to  the  implementation  of  boundary 
conditions.  Figure  2  shows  the  boundary  conditions  and  the 
mesh  that  was  used  in  the  calculations.  Once  the  nosing 
process  begins,  the  top  portion  of  the  tube  enters  the  die  zone. 
The  boundary  conditions  in  the  contact  region  are  stated  not 
only  by  prescribing  tractions  and  incremental  displacements, 
but  by  their  ratios.  The  radial  and  axial  velocities  of  the 
material  nodes  in  the  contact  region  are  not  independent,  and 
they  are  related  to  each  other  through  a  mathematical  ex¬ 
pression  for  the  geometrical  requirement  that  they  must  be 
actually  on  the  surface  of  the  die.  The  expression  is 

(r0+«)J+(Zo+*)1“('b  +  dr+<i)J+(Zo  +  Az+l>)J“/fJ. 

where  r0,  z0  are  radial  and  axial  positions  of  the  node  at  the 
present  configuration  at  time  f0,  and  dr,  dz  are  the  increments 
of  radial  and  axial  displacements.  R  is  the  radius  of  curvature 
of  the  nose  ogive.  See  Fig.  3. 

Rewriting  dr  and  dz  in  terms  of  velocities  u,  and  ut  as  dr  * 
ur dr  and  dz  -  utdr,  and  linearizing,  equation  (5)  results  in 
the  boundary  condition  given  by 

du,-^+j8  (6) 

a 

where 

^  (r0  +  a)  +wrdr  1 

(z+b)+n,d/  *  tan«’ 

_  [2(r0  +g)u, +  u?  +  dr  +  2(z0  +6)u.  -fi^d/l 
2l(r„+<i)+u,d/] 

Coulomb  friction  is  assumed  to  be  operative  between 
workpiece  and  die.  First,  we  prescribe  a  tangential  friction 
force  5  and  obtain  a  converged  solution.  Then,  we  compute 
the  generalized  nodal  forces.  With  angle  a  defined  in  equation 
(7),  we  then  are  able  to  compute  the  normal  component  N  and 
the  friction  coefficient  m  “  S/N  corresponding  to  the  initially 
prescribed  value  of  S.  If  the  computed  friction  coefficient  is 
not  what  is  intended,  we  then  modify  the  friction  force  S  and 
repeat  the  computational  procedure. 

Simulation  of  Cold  Nosing  of  105-mm  Shell  of  A1S1 
1018  Steel 

In  order  to  illustrate  the  application  of  the  FEM  analysis, 
simulations  were  made  for  cold  nosing  of  a  105-mm  shell. 
Some  cold  nosing  experiments  were  conducted  at  the  Battelle 
Columbus  Laboratories  (9),  and  the  dimensions  were  selected 
from  these  experiments,  as  shown  in  Table  1 .  The  specimens 
were  cut  from  AISI 1018  mild  steel  tubing,  and  its  flow  stress 
was  determined  at  Battelle  by  conducting  uniform  com¬ 
pression  tests.  The  actual  true  stress  versus  true  strain  curve 
(»  -  0  is  given  approximately  by  the  expression 

*-37.0(1 +50*)0*4  (8) 

where  *  is  the  flow  stress  at  an  effective  strain  I.  In  the 
computer  program,  we  have  constitutive  relationship  of  the 
form 

•-’’Mt)']  m 
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By  making  y  «  00,  and  Y  -  0,  the  relationship  given  by  (8)  is 
easily  incorporated. 

This  nonsteady-state  deformation  problem  was  analyzed  in 
a  step-by-step  manner  by  treating  it  quasi-linearly  during  each 
incremental  deformation.  Ten  uniform  elements  were  used  in 
the  axial  direction.  The  initial  step  size  for  each  incremental 
deformation  step  is  6  mm.  When  the  boundary  contact  region 
moves,  the  step  size  is  modified  in  such  a  way  that  a  free  nodal 
point  just  comes  in  contact  with  the  die  at  the  end  of  the  step. 

The  nosing  velocity  was  6.0  in. /min,  or  2.54  mm/s.  The 
uniform  velocity  field  of  2.54  mm/s  was  used  as  an  initial 
guess.  The  solution  obtained  from  the  previous  step  was  then 
used  as  an  initial  guess  for  the  subsequent  step.  The  com¬ 
putation  was  performed  for  each  step  solution  until  an  ac¬ 
curacy  of  the  error  norm  lAvl/M  s  0.00001  was  reached.  It 
took  8  to  10  iterations  to  reach  this  point  of  cor  ergence.  For 
each  step,  the  friction  forces  were  modified  three  to  four  times 
to  get  to  within  5  percent  of  error  of  the  desired  friction 
coefficient  value.  The  considered  offset  value  i0  described 
previously  was  0.0001 . 

The  calculated  load-displacement  curves,  with  friction 
coefficients  of  0,  0.05,  0.1,  0.15,  and  0.2,  respectively,  are 
shown  in  Fig.  4.  It  is  observed  that  the  experimentally 
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measured  curve  of  BetteiW  Laboratories  stress  well  with  the 
FEM  simulated  curves  for  a  -  0.1.  In  the  experiments,  the 
specimens  were  phoephated  and  coated  with  commercial  soap 
lubricant  before  nosing.  It  is  interesting  to  note  that  a  -  0.1  is 
typically  used  in  cold  forging  analysis  of  steel  specimens  with 
phosphate  coating  and  soap  lubrication. 

Figure  5  shows  the  elongation  in  shell  length  due  to  nosing 


at  various  friction  values.  At  low  friction  values,  we  have  a 
greater  increase  in  length  than  at  high  friction  values.  This 
predicted  trend  is  in  agreement  with  experimental  ob¬ 
servations  reported  in  the  literature.  As  seen  in  Fig.  S,  again 
the  experimentally  obtained  elongation  due  to  nosing  agrees 
very  well  with  finite-element  simulations  for  friction  con- 
dition  n  -  0.1. 

During  the  nosing  process,  the  shell  wall  thickness  in¬ 
variably  increases  with  increasing  penetration  in  the  die  and 
with  increasing  frictional  restraint  at  the  die-workpiece  in¬ 
terface.  Figure  6  illustrates  this  trend  and  gives  the  shell-wall 
thickness  distribution  as  predicted  by  the  FEM  simulation  at 
various  friction  conditions.  The  thickness  of  the  shdl-wall 
progressively  increased  over  the  length  of  the  shell  in  contact 
with  the  die,  with  the  increase  being  greatest  at  the  nose  tip.  It 
is  observed  in  a  small  region  near  the  nose  tip  that  the  in¬ 
creases  in  trail  thicknesses  are  about  the  same  for  various 
friction  values. 

For  different  frictional  conditions,  the  distributions  of 


Tmnanotiona  of  ttw  A8ME 


CONTACT  LENGTH  FROM  OIE  ENTRANCE,  AH* 

Fl0.11  Dta  pressura  dlstrlbulkm  In  th*  noatng  of  AIS1 1011  sMal  aMt 


axial  strain  t:,  circumferential  strain  «#,  thickness  strain  cr, 
and  effective  strain  l  are  plotted  in  Figs.  7  to  9  for  the  entire 
length  of  the  nosed  shell  for  various  values  of  coefficient  of 
friction.  The  sum  of  t,,  <„  and  t,  remains  zero  for  any 
amount  of  plastic  deformation,  and  the  incompressibility 
condition  is  satisfied. 

It  is  interesting  to  note  from  the  strain  distribution  plots 
that  axial  strain  «,  has  negative  values  near  the  main  body  of 
the  shell.  This  suggests  that  the  workpiece  is  shortened  near 
the  entrance  of  the  die.  The  lengthening  of  the  shell  takes 
place  in  the  remaining  two-thirds  of  the  nosed  portion. 

The  determination  of  the  proper  preform  profile  to  use,  in 
order  to  obtain  a  certain  finished  shape  after  nosing,  involves 
a  consideration  of  the  strains  existing  in  the  shell  nose.  Figure 
10  shows  progressively  how  the  strain  distributions  change 
during  the  nosing  process  for  frictionless  conditions.  .The 
general  trend  is  the  same  for  all  friction  conditions. 

The  die  pressure  distribution  was  calculated  from  finite- 
element  simulation  results  by  dividing  nodal  point  forces  over 
the  element  contact  area.  In  Fig.  11,  the  die  pressure 
distribution  is  plotted  along  the  die  contact  length  for  various 
friction  conditions  at  the  end  of  nosing.  The  yield  strength  of 
the  material  is  oy. 

It  can  be  teen  from  the  figure  that  friction  at  the  die- 
workpiece  interface  has  little  effect  on  die  pressure 
distribution.  It  is  also  noted  that  the  die  pressure  at  the  die 
entrance,  or  the  nose  base,  is  of  the  magnitude  of  the  yield 
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strength  of  the  material  at  the  nosing  limit.  This  indicates  that 
for  friction  formulations  in  FEM  simulation,  both  the 
Coulomb  Law,  F  ■  iiP,  and  the  Shear  Friction  Formulation, 
F  ■  mk  (k  is  shear  yield  stress)  will  give  similar  results.  In 
other  words,  the  values  of  a  and  m  in  shell  nosing  process  are 
comparable. 

The  limit  of  die  penetration  was  recognized  in  the  finite- 
element  simulation  by  the  fact  that  the  solution  diverges  after 
a  certain  amount  of  penetration.  The  limit  of  die  penetration 
depends  on  friction  along  the  interface,  and  the  effect  of 
friction  coefficients  on  this  limit  is  illustrated  in  Fig.  12.  As 
the  coefficient  of  friction  is  increased,  the  penetration  limit  is 
reduced,  according  to  the  computer  simulation. 

Investigation  of  the  Nosing  of  Shells  for  2024 
Aluminum 

In  order  to  make  more  extensive  comparisons  between 
theory  and  experiment,  laboratory  experiments  were  con¬ 
ducted  to  study  the  geometrical  changes  and  the  forces 
required  for  nosing  shells  of  aluminum  alloy  2024. 

The  geometry  of  the  die  shape  and  the  dimension  of  the  die 
and  the  workpiece  comprise  a  scaled-down  model  of  the  105- 
mm  shell.  The  specimen  dimensions  are  1.7  in.  (43.18  mm) 
outside  diameter,  3.6  in.  (91.44  mm)  length,  and  0.204  in. 
(5.18  mm)  wall  thickness.  Nosing  was  conducted,  using  wax 
as  a  lubricant,  with  a  speed  of  0.1  in./min.  Four  specimens 
were  tested,  but  pushed  to  various  depths  into  the  die.  The 
load  and  displacement  curves  were  recorded.  After  nosing, 
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the  specimens  were  cut  along  a  diametric  plane  and  elongation 
and  thickness  distribution  were  measured.  The  true  stress  and 
strain  relationship,  given  by  o  ■  10  (1  +  4.6*° 2n)  (ksi),  was 
used  in  the  flnite-demcnt  calculation.  The  over-all  visual 
comparison  of  the  deformed  shapes,  obtained  by  computation 
and  by  experiment,  is  shown  in  Fig.  13.  Good  agreement 
between  them  in  terms  of  total  height  and  wall  thickness 
change  is  evident.  Computer  simulations  were  made  for  two 
friction  conditions:  a  «  0.02  and  p  *  0.03.  The  calculated 
load-dispiacstnent  curves,  together  with  the  experimental 
result,  are  presented  in  Fig.  14.  It  is  observed  that  the  ex¬ 
perimental  remit  agrees  very  well  with  the  simulation  for 
friction  condition,  a  ”  0.02. 

The  elongation  of  specimens  due  to  nosing,  together  with 
the  elongation  predicted  by  the  finite-element  method,  are 
shown  in  Fig.  13.  Under  the  test  conditions,  the  largest  in¬ 
crease  in  length  was  about  2.3  percent.  Again,  the  agreement 

• 


19 


AXIAL  DISTANCE  PROM  BASE,  ma 

PI—  a  m  W I  —  _  m m  -—A  -I  m  MW  -  -  -  —  Si 

rig,  10  KMTMOPwMnT  Of  HiOioMft  ml  WW  uuvwlWI 


between,  the  experimental  results  with  finite-element 
simulations  for  friction  condition  p  -  0.02  are  good.  As  the 
depth  of  penetration  increased,  the  experimental  results 
moved  closer  towards  the  finite-element  predictions  for  the 
friction  condition  of  p  *  0.03.  This  is  because  the  speed  of 
nosing  was  slow  and  the  lubricant  was  squeezed  out  as  the 
specimen  penetrated  deeper  into  the  die.  Then  the  friction 
condition  between  shell  and  die  no  longer  remained  the  same 
as  at  the  beginning  of  nosing,  and  the  lengthening  was 
reduced.  It  is  also  noted  that  the  change  in  length  becomes 
more  sensitive  to  the  lubrication  condition  as  the  depth  of 
penetration  increased.  This  illustrates  the  importance  of 
lubrication  in  nosing. 

The  development  of  increase  in  wall  thickness  at  various 
depths  of  penetration  is  shown  in  Fig.  16.  The  finite-element 
predictions  for  friction  condition  p  *  0.02  are  also  presented 
for  comparison.  The  results  shown  in  Figs.  14,  13,  and  16 
demonstrate  that  agreement  between  theory  and  experiment  is 
excellent,  not  only  over-all  but  also  in  detail. 


The  finite-element  simulation  methods  were  developed  to 
study  nosing  of  shells  at  room  temperature.  In  the  finite- 
element  formulation,  the  nine-node  quadrilateral  dements 
with  quadratic  velocity  distribution  were  used.  In  order  to 
satisfy  the  geometrical  requirement,  the  boundary  conditions 
of  the  material  nodes  in  the  die  contact  region  were  im¬ 
plemented  through  a  mathematical  expression  which  relates 
the  axial  velocity  to  radial  velocity  during  the  nosing.  With  the 
advancement  of  the  die,  the  boundary  separating  the  contact 
region  from  the  unsupported  region  changes.  The  treatment 
of  moving  boundary  in  the  analysis  of  nosing  was  successfully 
implemented  in  the  finite-dement  program. 

Computer  simulations  were  performed  for  nosing  shells  of 
AISI  1018  steel  and  aluminum  alloy  2024.  The  solutions  were 
obtained  in  terms  of  load -displacement  curves,  thickness 
distribution,  elongation  and  strain  distributions.  The 
solutions  for  1018  sted  nosing  were  compared  with  ex¬ 
perimental  results  obtained  at  the  Battdle  Columbus 
Laboratories.  Experiments  of  aluminum  2024  nosing  were 
conducted  to  observe  metal  flow  characteristics  in  more 
detail. 

Comparisons  between  the  predictions  and  experimental 
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results  were  satisfactory,  not  only  over-all  but  also  in  detail 
for  both  materials. 
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An  Investigation  of  the  Shell 
Nosing  Process  by  the  Finite- 
Element  Method.  Part  2:  Nosing 
at  Elevated  Temperatures  (Hot 
Nosing) 

The  metal-forming  processes  of  shell  nosing  at  elevated  temperatures  were  analyzed 
by  the  finite-element  method.  The  strain-rate  effects  on  materials  properties  and  the 
flo  w  stress  dependence  on  temperatures  were  included  in  the  finite-element  analysis. 
A  thermodynamic  theory  of  visco-plasticity  based  on  rational  mechanics  was 
adapted  to  a  rigid-plastic  material  idealization.  An  implicit  scheme  is  used  for  the 
time  integration  of  heat  transfer  equations,  which  is  coupled  to  the  plasticity 
equations.  The  nine-node  quadrilateral  elements  with  quadratic  velocity 
distribution  were  used  for  the  workpiece,  and  four-node  quadrilateral  elements 
were  used  for  the  die  in  the  heat  transfer  analysis  and  temperature  calculations.  The 
coupled  analysis  of  heat  transfer  and  deformation  was  applied  to  the  forming  of 
AISI 1045  steel  shells.  Correlation  between  simulation  and  experimental  results  are 
good. 


Introduction 

Nosing  of  large  shells  (133  mm  and  above)  is  done  at  hot 
working  temperatures.  Prior  to  nosing,  the  end  portion  of  the 
tube  to  be  formed  is  heated  to  temperatures  of  between  1300 
to  1900*F  by  tubular  wound  induction  coils.  Then  the  shell  is 
nosed  by  forcing  it  into  a  suitably  shaped  die.  When  defor¬ 
mation  takes  place  at  high  temperature,  the  flow  stress  of  the 
material  is  a  function  of  temperature  and  strain  rate.  The 
temperatures  in  the  shell  wall  are  influenced  by  (1)  preheating 
conditions,  (2)  heat  generation  due  to  deformation  and 
friction,  and  (3)  heat  transfer  to  the  nosing  die.  Determination 
of  these  temperature  distributions  along  the  length  of  the  shell 
wall  is  important  and  necessary  in  order  to  predict  accurately 
the  flow  of  metal  and  the  nosing  load  in  the  hot  nosing. 

In  fact,  the  temperature  distribution  is  rather  critical  and 
much  effort  has  been  concentrated  on  the  preheating  con¬ 
ditions.  If  the  heating  is  gradual  and  slow,  the  unheated 
portion  of  the  tube  gains  heat  by  conduction.  Consequently, 
the  tube  tends  to  buckle  under  axial  load  during  nosing. 
Therefore,  preheating  is  done  at  a  rapid  rate  so  that  the  high 
temperature  does  not  travel  back  to  that  part  of  the  shell 
which  should  remain  relatively  cold. 

Until  recently,  most  of  the  temperature  calculations  during 
a  metal-forming  process  have  been  based  on  a  procedure  that 
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uncouples  the  problem  of  heat  transfer  from  the  metal 
deformation  problem.  The  approach  was  to  determine  the 
flow  characteristics  of  the  problem  either  experimentally  or  by 
calculations,  then  use  these  calculations  for  the  computation 
of  heat  generation  and  temperature  distribution. 

A  paper  by  Oden,  Bhandary,  Yagawa,  and  Chung  [1] 
presents  the  first  coupled  analysis  of  deformation  and  heat 
transfer  for  elasto-viscoplastic  materials  in  a  three- 
dimensional  rectangular  bar.  Rebelo  and  Kobayashi  [2]  have 
made  a  coupled  thermo-viscoplastic  analysis  in  solid  cylinder 
and  ring  compression.  In  a  previous  work  by  Lahoti, 
Subram anian,  and  Altan  [3),  the  heat-transfer  problem 
during  and  after  induction  heating  was  analyzed  by  the  Finite- 
difference  method.  We  present  here  the  finite-element 
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simulation  for  the  coupled  analysis  of  metal  deformation  and 
heat  transfer  during  hot  nosing.  The  results  are  compared 
with  the  experimental  work  by  Carlson  [4],  and  with  the 
confirmation  testing  of  hot  nosing  process  conducted  by 
Battellc  Columbus  Laboratories  [5].  The  effects  of  tern* 
perature  distribution  on  nosing  are  shown  and  discussed. 

Coupled  Thermo-  Viscopiastic  Analysis 

The  finite-element  equations  resulting  from  discretizing 
weak  forms  of.  the  equilibrium  equations  and  heat  balance 
equations  were  given  by  Rebelo  and  Kofcayashi  [2].  The  finite- 
element  method  for  the  deformation  analysis  is  described  in 
the  accompanying  paper  [6]. 

For  heat-transfer  analysis,  heat  balance  equations  upon 
finite-element  discretization  are  reduced  to  the  form 

Ct+/fT«Q,  (1) 


C  »  heat  capacity  matrix 
K  -  heat  conduction  matrix 
Q  *  heat  flux  vector 
T  -  vector  of  nodal  point  temperatures 
T  «  vector  of  nodal  point  temperature  rates 

The  heat  flux  vector  Q  in  equation  (1)  has  several  com¬ 
ponents,  and  is  expressed  with  interpolation  function  N  by 


)NdSc 


+  jS(  Aub(r, -  r„)N dS,  +  j  9/N dST  (2) 

The  first  component  is  the  contribution  of  the  net  heat 
generated  inside  the  deforming  body,  and  the  fraction  v  of  the 
plastic  work  rate  is  assumed  to  be  0.9.  The  second  component 
is  the  contribution  of  the  heat  radiation  from  the  en¬ 
vironment,  where  a  is  the  Boltzman  constant,  t  is  the 
radiation  coefficient,  and  T,  and  T,  are  environment  tem¬ 
perature  and  surface  temperature,  respectively.  The  third 
term  describes  the  heat  convection  between  body  surface  and 
environment  with  heat  convection  coefficient  h.  The  fourth 
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term  represents  the  contribution  of  the  heat  transferred 
between  workpiece  and  die  through  their  interface.  T„  and  Tw 
are  die  temperature  and  workpiece  temperature,  respectively, 
and  /iLub  is  the  heat  transfer  coefficient  from  the  lubricant. 
The  last  term  is  the  contribution  of  the  heat  generated  by 
friction  along  the  die-workpiece  interface,  qf  being  the 
surface  heat  generation  rate  due  to  friction. 

Coupling  of  deformation  analysis  and  temperature 
calculation  is  seen  in  the  scheme  of  time  integration  of 
equation  (1). 

The  scheme  requires  consistency  and  stability.  The  first 
requirement  is  satisfied  by  an  approximation 

T,+*  «T,  +  A/l(l  +  (3) 

where  /  and  Ar  are  time  and  its  increment,  respectively,  and  8 
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Table  I  Straa,  strain,  strain-rale,  temperature  relationships  of  AISI 1045  sted 

t£i  +  (  — )  J  (unitKsi) 

_____ 

Y  -  -7,2J88x  I0'*r*  +  4A912xlO~47*  -  0.938T  +  -64.6 

h  -  -0.1743 x  IO'*^  +  0.1026  x  10" 5  7*  -  0.19I7X  I0*J T  +  1.3 

y  -  -0.8682x10**^  +  0.5283 x  lO*4^  -  0.107  +  72.28 

At  <  «  0.2 

Y  -  -4.9439x  lO**?4  +  3.1693  x  10 “ 4 71*  -  0.697  +  313.36 

n  -  0.7859xl0-*r)  -  0.4496 x  10  *  5 7*  +  0.8633 x  I0*2r -  J.38 
y  -  -0.501  xKr'T13  +  0.3057  x  I0‘4r*  -  0.0627+  42.06 
At «  -  0.3 

Y  -  -6.3728x10**7'*  +  4.0«2x  10*47*  -  0.8831  T  +  648.61 

n  -  -  0.5858x10**  7"*  +  0.3396 x  10 "*7*  -  0.6397 x  10_2r  +  4.09 
y  -  - 0.8535 x  10 ‘‘r1  +  0.5 19810*  4  7*  -  0.10537'+  71.05 
At  t  m  0.4 

y  -  -  0.8535  x  10 ''r*  +  0.SI98  x  lO*^2  -  0.1053r+  71.05 

n  m  -0.2457x10**7*  +  0.1278 x  10* 5 7*  -  0.2046X  I0*2r  +  1.15 

y  -  -0.6617x10**7*  +  0.4094  x  10  *  5 7*  -  0.8441  xlO'2r+  5.79 

where  T is  temperature  in  *F.  Use  interpolation  to  get  flow  stress  for  any  f  <  0.4 
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is  a  parameter  varying  between  0  and  1.  For  unconditional 
stability  0  should  be  greater  than  0.5,  and  a  value  of  0.75  was 
chosen.  A  detail  on  an  implementation  of  the  algorithm  to 
include  the  coupling  of  the  solution  of  stress  equilibrium 
equations  with  heat  balance  in  the  finite-element  program  is 
given  elsewhere  [6]. 

Hot  Noslag  of  a  Small  Model  of  105-mni  Shell 

To  study  the  nature  of  the  nosing  process  in  the  forging  of 
shells,  Carlson  [4]  conducted  hot  nosing  on  a  small  mode)  of  a 
lOS-mm  shell.  Commercial  open-hearth,  cold-rolled  steel  was 
used  in  making  the  specimens.  Geometric  details  of  the 
specimens  are  given  in  Fig.  1 . 

The  die  that  was  employed  throughout  the  experiments  has 
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a  radius  of  curvature  of  11.85  in.  for  an  ogive  profile.  To 
simulate  the  amount  of  nosing  done  on  the  105-mm  shell,  1 .45 
in.  of  the  specimen  was  made  to  enter  the  die,  giving  a 
maximum  reduction  in  the  mean  diameter  of  about  33  per¬ 
cent. 

In  preparation  for  the  test,  the  die  was  heated  for  about  1  hr 
to  bring  it  to  a  temperature  of  1000’F.  The  testing  tem¬ 
perature  of  the  shell  was  from  1600*F  to  1900*F  at  the  tip.  To 
measure  the  temperature  distribution,  six  chromel-alumel 
thermocouples  were  placed  along  the  shell  model.  The 
measured  temperature  distributions  are  shown  in  Fig.  2. 

The  heating  of  the  shell  was  done  in  about  I  min  by  an 
induction  coil.  The  time  of  nosing  was  of  the  order  of  1  s. 

The  finite-element  simulations  were  carried  out,  using  the 
mesh  system  and  dimensions  shown  in  Fig.  3.  A  total  of  39 
nodes  and  6  nine-node  elements  were  used  for  the  workpiece. 
The  die  mesh  has  13  four-node  quadrilateral  elements  and  a 
total  of  52  nodes. 

The  material  properties  of  the  specimen  used  in  the  ex¬ 
periments  were  very  close  to  that  of  AISI  1045  steel.  Flow 
stress  expressions  of  1045  steel  were  found  from  experimental 
data  of  Allan  and  Boulger  [7],  We  rgderive  the  expressions 
and  reduce  them  to  fit  the  relationship,  as  shown  in  Table  1 . 

The  heat  transfer  characteristics  were  taken  from  standard 
handbooks: 

Kw  -  30.28  N/s  K  (workpiece  heat  conductivity) 

(pC)w  -  3.69  N/mm2  K  (p,  density;  C„,  workpiece  heat 
capacity) 

Kd  »  19  N/s  K  (die  heat  conductivity) 

(pC)  „  -  3.69  N/mm2  K  (p,  density;  Cd,  die  heat  capacity) 
h  -  0.0692  N/s  mm  K  (heat  convection  coefficient) 

“  32.88  N/s  mm  K  (heat  transfer  coefficient  from 
lubricant) 
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DISTANCE  FROM  NOSE  TIP,  mm 

Fig.  0  Towporoluro  distribution  oflor  noolng  witti  Inttlol  Up  torn- 
poraturo  ■  1200  K 

The  outer  boundaries  of  the  die  are  sufficiently  far  from  the 
interface,  and  a  constant  temperature  is  imposed  on  time. 

The  nosing  speeds  of  47  nun/s  (1.85  in./s)  and  41  mm/s 
(1.62  in./s)  were  simulated  by  the  finite-element  calculations. 
The  uniform  velocity  field  equivalent  to  the  nosing  speed  was 
used  as  the  initial  guess.  Step  size  was  2.5  mm. 

In  deformation  calculations,  an  average  of  nine  iterations 
was  needed  for  convergence,  plus  up  to  ten  more  to  get  the 
desired  friction  coefficient.  Normally  another  three  iterations 
were  needed  for  temperature  convergence. 

The  FEM  load-displacement  curves  for  hot  nosing  at 
various  friction  conditions,  with  initial  tip  temperature  at 
1825*F  (1269  K)  and  speed  1.62  in./s  (41  mm/s)  are  given  in 
Fig.  4.  It  can  be  seen  that  the  simulation  results  are  very  close 


Fig.  •  Comp  Arison  of  final  lampwrrtmo  distributions  lor  dtttofsnt 
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Fig.  10  Effect  of  temperature  on  maximum  ferae  during  norifiQ 
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to  those  of  the  experiments  for  friction  coefficient  p  *=  0.1.  In 
Fig.  5,  the  load  displacement  curves  are  obtained  by  FEM  for 
initial  tip  temperatures  at  1625*F  (1158  K),  1750‘F  (1227  K), 
and  1825*F  (1269  K),  respectively,  at  nosing  speed  of  M 
mm/s  and  n  «  0.1.  Figure  6  is  a  plot  simitar  to  Fig.  5,  with 
higher  nosing  speed  of  47  mm/s.  All  of  the  FEM  simulation 
results  showed  good  agreement  with  the  experimental  data  for 
a  -  0.1.  It  is  observed  that  the  nosing  load  increased  as  the 
workpiece  temperature  decreased. 

The  effect  of  temperature  on  shell  elongation  is  illustrated 
in  Fig.  7.  It  is  noticed  that  as  the  wall  temperature  increased, 
the  elongation  is  reduced  slightly. 

Figure  8  shows  the  temperature  distribution  before  and 
after  nosing,  for  initial  tip  temperature  of  1269  K  and  nosing 
speed  of  41  mm/s,  at  n  -  0.1  and  n  «  0.05.  At  the  end  of 
nosing,  the  FEM  predicted  a  temperature  drop  of  about  160  K 
on  the  nose  tip.  However,  at  the  nose  base,  a  temperature  rise 


4 


Transactions  of  tfto  ASME 


Fig.  12  Tiraperatura  distribution  of  186-mmpraloim  prior  to  noting 


Fig.  13  Mooli  system  tor  FEM  simulation  of  164mm  shsll  noolng 

of  about  50  K  was  expected  due  to  conduction  heat  from  the 
upper  portion  of  the  shell  to  the  lower  portion.  Friction 
condition  has  very  little  effect  on  the  temperature 
distribution.  This  means  that  the  heat  generated  due  to 
friction  is  small  compared  to  the  amount  generated  by 
convection  or  conduction.  With  lower  initial  temperatures  on 
the  workpiece,  the  results  in  temperature  distribution  after 
nosing  are  similar.  A  comparison  of  final  temperature 
distributions  for  different  nosing  speeds  is  shown  in  Fig.  9.  In 
the  region  near  the  nose  tip  where  large  deformation  oc- 


Fig.  14  Companion  of  load-disploesiiisnt  curvas  Irani  FBI 
calculation  and  experiment  lor  hot  noting  of  155mm  srioll 


curred,  the  final  temperatures  are  predicted  to  be  higher  with 
greater  nosing  speed.  In  the  region  near  the  nose  base,  the 
temperatures  are  the  same.  More  heat  is  generated  when 
nosing  speed  increases. 

Figure  10  shows  maximum  force  curves  during  nosing  and 
at  initial  tip  temperature.  The  considerable  scattering  of  the 
experimental  results  is  due  to  the  limitation  of  the  indicator 
diagram  in  depicting  maximum  force.  However,  a  consistent 
trend  was  observed  when  comparisons  were  made  between 
FEM  simulation  and  experimental  results.  Deviations  between 
computer  and  experimental  values  were  about  IS  percent.  The 
effect  of  nosing  speed  on  maximum  force  was  less  obvious 
from  FEM  simulation  than  from  experiments.  This  may  be 
due  to  the  fact  that  the  rate-sensitivity  constant  used  in  the 
simulation  is  different  from  that  of  the  material  used  in  ex¬ 
periments. 

Hot  Nosing  of  155>nun  Shell 

Confirmation  testing  of  the  hot  nosing  process  was  con¬ 
ducted  by  Batteile  Laboratories  at  Chamberlain  Manufac¬ 
turing  Corp.  in  order  to  study  the  material  flow  of  a  selected 
preform  at  hot  nosing  temperatures.  The  projectile  was  that 
of  s  ’5-mm  shell,  and  the  material  was  modified  A1SI  1045. 

Ptctorm  dimensions  were  measured  and  recorded.  Pre¬ 
form  temperatures  were  measured,  using  imbedded  ther¬ 
mocouples  at  six  places  in  the  region  near  the  nose  tip.  Details 
of  the  preform  design  are  given  in  Fig.  11.  The  preform  was 
heated  to  1400*F  at  the  tip,  using  the  induction  coil.  The 
temperature  distribution  prior  to  nosing  is  shown  in  Fig.  12. 
The  nose  die  has  10.5-in  OD  and  is  enclosed  by  a  water-jacket; 
typically,  the  water  in  the  jacket  is  nearly  boiling.  The  die 
temperature  is  estimated  at  about  400*F. 

Water-based  graphite  was  used  as  the  lubricant  in  the  test. 
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The  ram  speed  of  the  800-t  mechanical  press  was  2.3  in.  (63.3 
mm)  per  second  in  carrying  out  the  nosing  process. 

The  mesh  system  for  the  finite-element  simulation  was  built 
with  39  nodes  and  6  nine-node  elements  in  the  workpiece,  and 
36  nodes  and  39  four-node  elements  in  the  die.  The  initial  grid 
is  shown  in  Fig.  13.  Total  penetration  of  the  shell  into  the 
nosing  die  was  about  223  mm,  obtained  in  14- mm  steps.  A 
constant  temperature  was  imposed  on  the  outer  boundaries  of 
the  die. 

Flow  stress  expressions  for  AIS1  1043  and  the  heat  transfer 
characteristics  were  the  same  as  described  in  the  previous 
section. 

The  load-displacement  curves  that  were  obtained  from 
finite-element  simulations  for  various  friction  coefficients  are 
compared  with  experiments  in  Fig.  14.  The  ring-compression 
test  gave  a  value  of  p  -  0.09,  and  from  FEM  results  this  value 
was  considered  to  be  too  high.  A  value  of  p  ■  0.02  or  smaller 
gives  simulation  results  closer  to  that  in  the  experiments.  Near 
the  end  of  nosing,  the  experiments  showed  a  steeper  increase 
in  the  load  than  did  the  FEM,  possibly  because  the  lubricant 
was  squeezed  out,  thereby  increasing  friction. 

The  temperature  distributions  of  the  shell  outer  surface  at 
different  penetration  depths  are  presented  in  Fig.  13.  As  the 
contact  between  die  and  workpiece  increases,  the  temperature 
of  the  workpiece  in  the  contact  region  gradually  decreases  due 
to  heat  flow  from  the  workpiece  to  the  colder  die.  The  test 
results  showed  product  temperatures  of  between  900  and  1000 
K,  which  agree  with  the  average  temperature  distribution 
from  the  FEM. 


DISTANCE  T  AON  NOSE  TIA.  mm 

Pig.  18  Strain  distributions  at  penetration  ol  220  mm 


Pig.  19  Comparison  of  deformed  shapes  between  FEM  and  ex¬ 
periment 

Temperature  rises  as  expected  around  the  contact  areas  in 
the  die,  and  the  differences  between  workpiece  and  die 
temperatures  decrease  substantially  during  the  course  of 
deformation. 

A  comparison  of  thickness  distribution  along  the  shell  is 
given  in  Fig.  16.  Thickness  distributions  at  several  penetration 
depths  are  predicted  by  the  FEM  model.  An  increase  of 
penetration  depth  increases  the  thickness  for  a  given  point. 
Also,  the  region  of  wall  thickening  extended  towards  the  nose 
base.  A  consistent  trend  was  observed  when  comparisons  were 
made  between  computer-calculated  wall  thickness  and  ex¬ 
perimental  results. 

The  133-mm  shell  length  increases  after  hot  nosing  under 
good  lubrication,  as  shown  both  by  experiments  and  by  FEM 
simulation.  Figure  17  gives  shell  elongation  versus  die 


6 


Tranaaetlona  of  tho  ASME 


I 


penetration  during  hot  nosing.  An  increase  in  length  of  8 
percent  is  expected  after  penetration  of  220  mm.  Recall  that 
only  4  percent  elongation  was  expected  at  the  end  of  nosing 
for  a  105-mm  shell.  This  indicates  more  elongation  for  larger 
reduction  in  mean  diameter  in  shell  nosing. 

In  Fig.  18  we  plotted  the  distributions  for  axial  strain  ez, 
radial  strain  tr,  and  circumferential  strain  u  st  die 
penetration  of  220  mm.  The  compatibility  is  satisfied  as  the 
sum  of  remains  zero.  The  axial  strain  e.  does  not 

show  negative  value  near  the  nose  base  as  it  did  in  the  cold 
nosing  simulation.  This  suggests  that  lengthening  takes  place 
in  the  entire  nosed  portion  of  the  shell. 

The  deformed  shape  of  the  nose  from  FEM  simulation  is 
compared  to  that  of  the  experiment  in  Fig.  19.  Since  the 
preforms  are  not  exactly  the  same  in  the  region  near  the  nose 
tip,  the  FEM  predictions  did  not  match  the  experiments. 
Additionally,  the  FEM  model  showed  a  nosing  limit  beyond 
which  material  can  no  longer  support  the  nosing  load,  and  a 
change  of  deformation  mode  (e.g.,  buckling)  may  occur. 
Therefore,  the  maximum  depth  of  penetration  was  also 
different  for  the  FEM  model.  The  discrepancy  may  arise  from 
the  fact  that  the  flow  stress  expressions  implemented  in  FEM 
calculations  could  be  different  from  those  of  the  materials 
used  in  experiments. 


Summary  and  Conclusions 

In  this  investigation  the  finite-element  simulation  methods 
were  developed  to  study  nosing  of  shells  at  elevated  tem¬ 
peratures.  The  deformation  analysis  and  heat  transfer 
analysis,  including  dies,  were  coupled  together  for  an  iterative 
simultaneous  solution.  The  nine-node  quadrilateral  elements 
with  quadratic  velocity  distribution  were,  used  for  the 
workpiece,  and  four-node  quadrilateral  elements  were  used  in 
the  die  for  the  heat  transfer  analysis  and  temperature 
calculations. 

The  equations  that  describe  the  deformation  mechanics  are 
nonlinear  and  were  solved  by  Newton -Raphson’s  iterative 
technique.  An  implicit  scheme  was  used  in  the  time  in¬ 
tegration  of  heat  transfer  equations  which  provides  good 
accuracy  for  a  reasonably  large  range  of  time  steps.  The  die 
part  was  remeshed  after  each  deformation  step. 

Finite-element  programs  were  written  in  FORTRAN  IV  for 
a  CDC  7600  computer.  Listing  of  these  programs,  including 
instructions  for  input  data  preparation,  are  given  in  reference 
(81- 

It  is  important  to  determine  the  precise  shapes  of  the 
profiles  of  shells  (thickness  distribution)  in  the  original 
cylindrical  forgings  to  produce  the  prescribed  shapes  of  shells 
in  their  finished  co*  dition.  For  this  requirement,  the  deter¬ 


mination  of  deformation  characteristics,  in  particular  the 
thickness  changes,  as  functions  of  process  variables  is 
essential.  Some  factors  which  influence  deformation 
characteristics  significantly  are  the  temperature  distribution 
along  the  cylinder,  friction  along  the  die,  and  pressing  speeds. 

The  coupled  analysis  of  heat  transfer  and  deformation  was 
performed  by  simulating  nosing  of  a  155-mm  shell  of  AISI 
1045  steel.  Experimental  results  were  available  from  Battelle 
Laboratories,  and  correlation  was  good. 

By  using  the  finite-element  programs  developed  in  this 
research,  we  are  able  to  obtain  :uch  detailed  deformation 
characteristics  as:  load-displacement  curves,  thickness 
distribution,  elongation,  strain  and  strain-rate  distributions, 
and  temperature  distributions.  They  are  valuable  for  the 
control  of  the  nosing  process,  design  of  equipment,  and 
preparation  of  workpieces. 
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ABSTRACT 


The  finite-element  method  for  rigid-plastic  materials  is  used  for  the 
analysis  of  nosing  of  aluminum  alloy  2024  tubular  preforms  with  uniform  wall 
thickness  at  room  temperature.  The  detailed  mechanics  were  obtained  by 
computer  simulation  and  some  of  the  results,  such  as  load  displacement 
curves,  elongation  and  thickness  distributions,  were  compared  with  experi¬ 
mental  results.  Then,  assuming  that  the  strain  distributions  computed  by 
the  finite-elemdnt  method  are  valid  for  preforms  having  nonuniform  thick¬ 
nesses  initially,  two  types  of  preform  shapes  were  desiqned  to  obtain  uniform 
wall  thickness  after  nosing.  Experimental  results  of  nosing  these  preforms 
showed  that  the  approach  to  the  problem  is  satisfactory. 
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ABSTRACT 

A  new  approach  to  preform  design  is  introduced,  as  a  unique  application 
of  the  finite  element  method  to  the  problems  in  metal  forming.  The  concept 
involved  in  the  approach  is  to  trace  backward  the  loading  path  in  the  actual 
forming  process  from  a  given  final  configuration  using  the  finite  element 
method.  The  concept  was  tested  with  the  problem  of  simple  compression  of  a 
cylinder  and  the  method  was  applied  to  preform  design  in  shell  nosing. 
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1 ♦  INTRODUCTION 


One  of  the  most  important  aspects  of  the  closed  die  forging  process  is 
the  design  of  preforms.  An  example  of  preform  design  for  steel  finish  forg¬ 
ings  of  various  H-shapes  is  shown  in  Fig.  la.  Preform  design  involves  the 
determination  of  number  of  preforms  and  design  of  the  shapes  and  dimensions 
of  each  preform.  According  to  Biswas  and  Knight  [1],  the  bulk  of  the  work  to 
determine  the  preform  shapes  for  complex  forging  components  has  been  carried 
out  in  Germany  and  the  Soviet  Union,  as  reported  by  Spies  [2].  Chamouard  [3] 
in  France,  and  Akgerman  et  al.  [4]  in  the  U.S.  dealt  with  some  practical  pro¬ 
blems  of  preform  design  for  certain  forging  cross-sections.  The  present  prac¬ 
tice  to  this  complex  problem  is  to  computerize  the  design  calculations  required 
for  the  designer's  decisions,  based  on  qualitative  guidelines  which  have  been 
derived  mainly  from  experience  or  experimental  studies. 

Preform  design  problem  is  also  involved  in  shell  nosing.  Shell  nosing  is 
the  process  of  forming  an  ogive  nose  at  an  open  end  of  a  tubular  part  (preform) 
by  pressing  the  tube  into  a  contoured  die.  Wall  thickness  increases,  but  the 
shell  may  elongate  or  shrink  in  length  during  nosing.  Because  shell  specifi¬ 
cations  usually  require  certain  wall  thickness  distributions  after  nosing,  the 
problem  is  to  design  the  preform  with  which  the  desired  shape  after  the  nosing 
operation  can  be  achieved.  An  example  of  such  a  preform  is  shown  in  Fig.  lb. 

The  approach  to  the  preform  design  in  shell  nosing  has  been  more  quanti¬ 
tative,  because  the  number  of  variables  involved  in  the  determination  of  pre¬ 
form  shapes  in  shell  nosing  is  fewer  than  those  for  preform  shapes  in  general 
forgings.  Nadai  [5]  outlined  an  approach  for  designing  the  preforms  for 
shells  with  uniform  wall  thickness  after  nosing,  based  on  the  strains  in  the 
shell  nose.  Carlson  [6]  extended  Nadai 's  approach  and  suggested  the  method 
for  determining  the  original  shell  profile.  Based  on  Carlson's  methqd,  Lahoti 
et  al .  [7]  developed  a  computer  program  for  the  preform  design. 
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In  this  paper  a  new  approach  to  preform  design  is  introduced.  The  concept 
involved  in  the  approach  is  to  trace  backward  the  loading  path  in  the  actual 
forming  process  from  a  given  final  configuration  by  the  finite  element  method. 
The  concept  was  tested  with  the  problem  of  simple  cylinder  compression  and  the 
method  was  applied  to  preform  design  in  shell  nosing.  Because  the  computa¬ 
tional  technique  uses  the  finite  element  method,  a  brief  description  of  the 
finite  element  method  is  given,  and  the  new  approach  to  the  problem  follows. 

2.  FINITE-ELEMENT  METHOD 

The  concept  of  the  finite-element  method  is  one  of  discretization.  Ttfe 
finite-element  model  is  constructed  in  the  following  manner.  A  number  of 
finite  points  is  identified  in  the  domain  of  the  function,  and  the  values  of 
the  function  and  its  derivatives,  when  appropriate,  at  these  points  are  speci- 
fied.  These  points  are  called  nodal  points.  The  domain  of  the  function  is. 
represented  approximately  by  a  finite  collection  of  subdomains  called  finite 
elements.  The  domain  then  is  an  assemblage  of  elements  connected  together 
appropriately  on  their  boundaries.  The  function  is  approximated  locally  with¬ 
in  each  element  by  continuous  functions  which  are  uniquely  described  in  terms 
of  the  nodal  point  values  associated  with  the  particular  element.  The  path  to 
the  solution  of  a  finite-element  problem  consists  of  five  specific  steps: 

(a)  the  problem,  (b)  the  element,  (c)  the  element  equation,  (d)  the  assemblage 
of  element  equations,  and  (e)  the  numerical  solution  of  the  global  equations. 

The  basis  of  finite-element  metal  flow  simulation,  using  the  variational 
approach  is  to  formulate  proper  functionals,  depending  upon  specific  constitu¬ 
tive  relations.  The  solution  of  the  original  boundary  value  problem  is  ob¬ 
tained  by  the  solution  of  the  dual  variational  problem  where  the  first-order 

variation  of  the  functional  vanishes.  Choosing  an  approximate  interpolation 
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function  for  the  field  variable  in  the  elements,  the  functional  is  expressed 


i 


locally  within  each  element  in  terms  of  the  nodal  point  values.  The  local- 
element  equations  are  then  assembled  into  the  overall  problem.  Thus,  the 
functional  is  approximated  by  a  function  of  global  nodal  point  values.  The 
condition  for  this  functional  to  be  stationary  results  in  the  stiffness  equa¬ 
tions.  For  metal-forming  problems  the  stiffness  equations  are  nonlinear  and 
solved  by  the  Newton-Raphson  iterative  method. 

For  rigid-plastic  materials,  obeying  the  Mises  yield  criterion  and  its 
associated  flow  rule,  the  incompressibility  constraint  can  be  removed  by  intro¬ 
ducing  either  the  Lagrange  multiplier  or  a  penalty  function  into  the  functional. 
If  the  body  contains  a  rigid  region,  the  stresses  in  the  region  cannot  be 
determined.  This  difficulty  is  handled  by  considering  an  offset  of  the  effec¬ 
tive  strain  rate,  which  is  several  orders  of  magnitude  smaller  than  the  aver¬ 
age  strain  rate  in  the  deforming  zone.  The  deviatoric  stresses  ?re  assumed  to 
vary  linearly  from  zero  to  the  flow  stress  if  the  effective  strain  rate  is 
smaller  than  this  offset  value.  This  stress-strain  rate  relation  in  the  "nearly 
rigid"  region  modifies  the  functional  accordingly. 

For  rigid-viscoplastic  materials,  the  constitutive  equations  derived  from 
the  viscoplastic  theory  are  similar  to  those  of  rigid-plastic  materials  except 
that  the  yield  stress  is  a  function  of  strain  and  strain-rate.  Applications 
of  these  methods  to  various  metal  forming  problems  are  given  in  a  review  on 
the  finite  element  method  and  metal  forming  process  modeling  [8],  Recently, 

Oh  [9],  at  the  Battel le  Columbus  Laboratories  has  accomplished  a  significant 
development  of  a  user-oriented  general  purpose  FEM  program,  namely  a  rigid- 
viscoplastic  FEM  Code  ALPID.  The  special  features  of  the  program  include: 

(1)  use  of  higher  order  elements,  (2)  general  description  of  die  and  automatic 
handling  of  die  boundary  conditions,  and  (3)  automatic  generation  of  Initial 
guess  for  Iteration  procedure.  The  program  ALPID  was  used  for  the  computation 
throughout  the  present  investigation. 
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3.  METHOD  OF  APPROACH 

A  method  of  approach  is  illustrated  in  Fig.  2.  At  time  t  *  tQ,  the  geo¬ 
metrical  configuration  Xg  of  a  deforming  body  is  represented  by  a  point  Q. 

The  point  Q  is  arrived  at  from  the  point  P  whose  configuration  is  given  as 
Xq_i  at  t  s  tg_i ,  through  the  displacement  field  Ug_^  during  a  time  step  At,  namely, 
xg  *  Xg_.j  +  Uq_i  .  Therefore,  the  problem  is  to  determine  Ug_.j,  based  on  the 
information  (xg)  at  point  Q.  The  solution  scheme  is  as  follows:  taking  a 
loading  solution  Ug  (forward)  at  Q,  a  first  estimate  of  P  can  be  made  accord¬ 
ing  to  p^  *  Xg  -  Ug.  Then,  the  loading  solution  u^j  can  be  calculated  on 
the  basis  of  the  configuration  of  p^  with  which  the  configuration  xQ  at  Q 
can  be  compared  with  p^  +  Ug^j  *  Q^.  If  Q  and  are  not  sufficiently 
close  to  each  other,  then  p^  can  be  estimated  by  p^  =  xQ  -  u^j.  The  sol¬ 
ution  for  loading  at  p^  is  then  Ug^j  and  the  second  estimate  of  the  con-  . 
figuration  Qv  '  -  p'  '  +  Ug  (  can  be  made.  The  iteration  is  carried  out  until 
q(n)  _  p(n)  +  ^n)  j,ecomes  sufficiently  close  to  Q.  Displacement  solutions 
are  computed  by  the  finite  element  method.  The  sufficient  condition  for  con¬ 
vergence  is  that  the  body  is  plastically  deforming  at  all  times,  if  the  bound¬ 
ary  condition  does  not  change  at  P  and  Q.  We  call  this  calculation  scheme 
"Backward  Tracing  by  FEM."  So  far,  the  concept  includes  only  the  geometry  of 
deforming  materials,  and  some  questions  as  to  its  applicability  to  real  pro¬ 
blems  arise  immediately.  Some  of  them  are:  (i)  how  to  include  the  effect  of 
material  properties,  such  as  work-hardening,  strain-rate,  and  temperature 
effects;  (ii)  how  to  modify  the  method  when  the  boundary  conditions  change; 
and  (ill)  the  problem  of  solution  uniqueness. 

Although  complete  and  satisfactory  answers  to  these  questions  await  fur¬ 
ther  investigations,  a  test  with  axisymmetric  compression  of  a  cylinder  demon¬ 
strates  that  the  concept  does  work.  * 
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A  circular  cylinder  of  height-to-diameter  ratio  equal  to  1  is  compressed 
between  two  parallel  platens  with  friction  (friction  factor,  m  =  0.5)  to  10 
percent  reduction  in  height.  The  material  is  assumed  to  be  rigid-perfectly 
plastic.  In  Fig.  3  the  workpiece  configuration  at  10  percent  reduction  is 
shown  by  dotted  lines  as  the  result  of  rigid-plastic  finite-element  simula¬ 
tion  of  loading.  Now,  the  backward  tracing  scheme  was  applied  to  the  known 
configuration  at  10  percent  reduction  in  height.  After  the  10  percent  axial 
displacement  was  traced  backward,  the  result  indeed  gave  the  circular  cylinder 
as  the  preform. 

For  both  forward  loading  and  backward  tracing,  a  step  size  of  1%  reduc¬ 
tion  in  height  was  used.  Because  of  the  use  of  Newton-Raphson  method  in  the 
finite-element  program,  iteration  is  involved  in  obtaining  a  loading  solution 
(forward  iteration),  and  another  iteration  is  necessary  in  the  backward  trac¬ 
ing  scheme  (backward  iteration).  Numbers  of  iterations  for  each  step  of  de¬ 
formation  were  6-7  and  about  3  for  forward  and  backward  procedures  respec¬ 
tively.  Closeness  of  the  two  configurations  and  Q  was  measured  by  |xQ  -  x£ 
and  the  limit  was  set  as  0.005  (mm),  where  Xq  represents  the  coordinates  of 
the  configuration  of  Q^. 

It  should  be  noted  that  the  process  reversal  (reverse  loading)  by  revers¬ 
ing  the  direction  (sign)  of  the  boundary  conditions  does  not  trace  the  loading 
path  backward,  because  the  process  reversal  is  a  different  boundary  value  pro¬ 
blem  from  that  of  loading.  Figure  3  shows  that  clearly  the  shape  that  is 
obtained  by  reverse  loading  did  not  result  in  the  original  circular  cylinder. 

4.  APPLICATION  TO  SHELL  NOSING 

Preform  design  in  forging  is  complex  and  generally  involves  many  problems 

in  applying  the  new  approach.  Less  complex  is  the  preform  design  in  shell 

% 

nosing,  because  the  number  of  parameters  in  defining  preform  shapes  is  fewer. 
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In  shell  nosing,  two  simple  preform  configurations  can  be  realized  for  a  uni¬ 
form  wall  thickness  after  nosing.  The  wall  thickness  distributions  required 
in  the  preform  can  be  accommodated  by  varying  either  the  inner  or  the  outer 
diameters  of  the  preform.  It  is  not  difficult  to  see  that  different  preform 
configurations  can  be  arrived  at  by  imposing  different  sequences  of  boundary 
conditions  (such  as  where  and  when  the  nodal  point  in  contact  with  the  die 
should  be  freed  from  the  die)  during  the  backward  process.  This  aspect  of  the 
problem  is  first  investigated  for  non-workhardening  materials.  Then,  using 

work-hardening  materials,  the  procedure  to  include  the  effect  of  work-hardening 

».  ■ 

is  introduced  in  the  new  approach. 

Rigid-perfectly  plastic  materials: 

With  reference  to  Fig.  4,  specifications  of  the  final  nosed  shell  con¬ 
figuration  are  given  by  (unit;  mm) 

Lf  *  48.84  '  rQ  *  21.526  bf  *  14.47 

(Lq  *  47.26)  r^  *  14.648  hf  *  6.SB  (uniform) 

R  =  259  (aQ  =  237,  bQ  =  18.75) 

Coefficient  of  friction  at  the  die-workpiece  interface  is  assumed  to  be  u  *  0.05. 
For  the  finite  element  calculation,  40  four  node  elements  and  63  nodal  points 
are  used;  The  criterion  for  controlling  the  boundary  condition  during  the 
backward  tracing  depends  on  the  preform  shape  under  consideration. 

For  the  preform  shape  with  constant  outer  diamter  (Type-0  preform), 

r11  *  r  for  all  n,  where  r11  is  the  radius  of  a  nodal  point  located  along  the 
0  0  0 

outer  surface  of  the  shell.  The  boundary  condition  for  this  type  of  preform 
is  controlled  in  such  a  way  that  nodes  which  are  in  contact  with  the  die 
(ro  <  ro>  are  free^  from  the  die  when  condition  that  r£  3  rQ  is  reached 
during  the  backward  tracing. 

For  the  preform  shape  with  constant  inner  diameter  (Type-I  preform). 
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the  criterion  for  changing  the  boundary  condition  is  that  as  soon  as  rP  =  r.. 
for  the  nodes  r^  <  r^,  the  boundary  condition  of  the  corresponding  outside 
node  is  changed  from  that  of  die-contact  to  the  force  free  condition,  where 
rP  is  the  radius  of  a  nodal  point  along  the  inner  surface  of  the  shell. 

The  above  criteria,  however,  result  in  some  errors  because  the  deforma¬ 
tion  of  the  shell  is  not  confined  within  the  die,  but  extends  beyond  the  die- 
entrance.  In  order  to  reduce  the  error,  the  criteria  for  change  of  the 
boundary  condition  are  modified  as 

r£  *  rQ  -  AQ  for  Type-0  preform 

r?  8  r.  -  for  Type- 1  preform 

The  amount  A  in  the  modified  criteria  can  be  estimated,  for  example,  from  the 
loading  simulation  of  nosing  with  a  preform  of  uniform  wall  thickness. 

Using  the  modified  criteria  for  controlling  the  boundary  condition,  the 
computation  of  backward  tracing  was  carried  out  taking  2  mm  penetration  (in 
loading)  as  one  step  (total  22  steps  to  complete  the  calculation).  Solution 
convergence  needed  about  6-7  iterations  and  for  completing  one  step  backward 
tracing  1  -  2  iterations  for  the  case  where  boundary  condition  does  not  change, 
2-3  iterations  are  required  for  the  case  where  change  in  the  boundary  condi¬ 
tion  occurs.  Convergence  for  backward  tracing  was  determined  by  closeness  of 
the  solution  within  the  limit  of  0.005  (mm).  The  two  types  of  preforms 

determined  by  the  method  are  compared  with  the  final  nose  configuration  in 
Fig.  5.  In  Fig.  5,  the  lengths  C®  and  for  Type-0  and  Type-I  preforms, 
respectively,  indicate  the  portion  of  the  preform  which  becomes  the  nosed 
shell  length  after  nosing.  Because  the  shell  elongates  during  nosing,  the 
preforms  are  shorter  than  the  nosed  shell.  It  can  also  be  observed  that  the 
Type-0  preform  is  slightly  shorter  than  the  Type-I  preform.  This  is,because 
the  amount  of  deformation  Involved  in  nosing  is  more  for  Type-0  than  for 
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Type- I,  and  the  total  elongation  is  larger  for  the  Type-0  preform.  The  thick¬ 
ness  distributions  in  the  two  preforms  are  compared  in  Fig.  6.  The  distribu¬ 
tions  are  almost  the  same  for  both  preforms,  except  that  the  thicknesses  tn 
Type-0  are  slightly  smaller.  This  is  reasonable  based  on  the  amount  of 
deformation  involved,  and  from  the  condition  of  volume  constancy. 

Work-hardening  materials: 

The  following  specifications  of  the  nosed  shall  are  used: 

material  =  A1  alloy  2024,  annealed,  with  the  stress- 
strain  curve  given  by 
5  *  68.9  (1.0  +  4.6(c)0,29)  Mn/m2 

Die  geometry  (mm):  same  (R  *  259,  a^  »  237,  b  *  18.77) 

o  o 

Shell  dimensions  (mm):  Lf  *  36.84,  bf  *  16.369, 

rQ  =  25.526,  r^  =  14.648,  hf  *  6.88 

Die-workpiece  interface  conditions:  u  s  0. 

The  shell  configuration  (in  terms  of  total  penetration)  is  different  from  the 
previous  non-workhardening  case.  The  reason  for  the  use  of  different  nosed 
shell  specifications  is  that  there  exists  a  limit  beyond  which  nosing  cannot 
be  continued  and  cylindrical  part  of  the  shell  begins  to  deform.  This  limit 
depends  on  the  work-hardening  characteristics  of  the  material,  and  the  present 
specifications  are  selected  within  this  limit. 

In  order  to  include  the  work-hardening  effect  in  preform  design,  it  is 
necessary  to  know  the  strain  distributions  in  the  nosed  shell,  in  addition  to 
the  geometrical  configurations,  but  the  strain  distributions  in  the  nosed 
shell  depend  on  the  preform  shapes,  and  are  not  known  (or  cannot  be  specified). 
In  Fig.  7  the  procedures  to  take  into  account  the  work-hardening  effect  in  the 
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present  approach  are  illustrated. 

Applying  the  backward  tracing  scheme  for  non-work-hardening  materials 
explained  in  the  previous  section,  the  two  types  of  preforms  can  be  determined. 
With  these  preforms  the  loading  solutions  are  obtained  for  a  work-hardening 
material.  The  results  provide  the  strain  (effective  strain)  distributions 
and  the  nosed  shell  configurations.  The  configurations  are  very  close  but  not 
identical  to  the  specifications.  These  configurations,  mainly,  the  thickness 
distributions,  are,  then,  corrected  exactly  according  to  the  specifications, 
maintaining  the  same  effective  strain  distributions.  With  these  strain  dis¬ 
tributions,  the  application  of  the  backward  tracing  scheme  results  in  the 
preform  shapes  for  a  work-hardening  material.  At  each  step,  the  strain  and 
corresponding  stress  are  also  traced  backward  following  the  stress-strain 
curve  of  the  material.  It  should  be  mentioned  that  the  first  computational 
procedure  in  Fig.  7  may  be  eliminated  by  using  approximate  solutions  [10]  for 
preform  shapes  of  non-work -hardening  materials. 

The  number  of  iterations  required  was  about  the  same  as  for  non-work¬ 
hardening  materials.  However,  solution  convergence  is  somewhat  sensitive  to 
the  size  of  time  step,  and  a  reduction  of  step  size  to  1  mm  penetration  was 
necessary.  This  step  size  required  total  30  steps  for  completing  the  computa¬ 
tions. 

The  effective  strain  distributions  in  the  nosed  shell  obtained  from  load¬ 
ing  solutions  for  both  types  of  preforms  are  shown  in  Fig.  8.  The  strains  are 
larger  for  Type-0  preform  than  for  Type-I  preform.  The  preforms  for  the  work¬ 
hardening  material  were  determined  by  applying  the  backward  tracing  scheme  to 
the  nosed  shell  with  the  strain  distributions  shown  in  Fig.  8,  and  the  speci¬ 
fied  geometrical  dimensions.  The  effective  strains  in  the  preforms  were  in¬ 
deed  negligible  (within  0.1%  strain),  when  the  backward  tracing  calculation 
was  completed.  The  thickness  variations  in  the  preform  of  the  work- 
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hardening  material  is  compared  with  that  of  non-work -hardening  materials  in 
Fig.  9(a)  for  Type-0  preform  and  in  Fig.  9(b)  for  Type-I  preform.  Fig.  9 
shows  that  the  preform  shapes  for  work-hardening  and  non-work -hardening  mater¬ 
ials  are  almost  identical  to  each  other  for  both  preform  types.  It  may  be  due 
to  the  fact  that  the  total  penetration  of  the  shell  into  the  die  is  so  small 
that  the  work-hardening  effect  is  negligible  for  this  specific  case.  Never¬ 
theless,  it  suffices  to  show  the  procedure  to  include  the  work-hardening 
characteristics  of  the  material  in  the  backward  tracing  technique. 

5.  CONCLUDING  REMARKS 

A  new  technique,  called  backward  tracing  scheme,  was  devised  for  preform 
design.  This  is  a  unique  application  of  the  finite  element  method  to  the 
problems  in  metal  forming.  It  was  demonstrated  that  the  technique  can  be 
applied  to  preform  design  in  shell  nosing.  It  was  also  revealed  that  one  of 
the  critical  aspects  for  further  development  of  the  technique  concerns  with 
the  treatment  of  the  boundary  conditions.  The  boundary  conditions  for  non¬ 
steady  state  forming  processes  change  during  loaning,  and  the  way  in  which 
this  change  occurs  depends  on  the  preform  shapes.  Thus,  the  criterion  for 
the  selection  of  desired  (or  optimum)  preform  are  needed,  and  then  the  con¬ 
dition  for  controlling  the  boundary  conditions,  corresponding  to  a  selected 
preform  must  be  developed. 

A  challenge  to  this  new  approach  for  further  development,  thus,  lies  in  its 
application  to  preform  design  in  forging. 
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ABSTRACT 

Two  approximate  solutions  for  preform  design  in  shell  nosing  are  derived. 

The  degree  of  approximation  was  examined  by  the  accurate  numerical  solution 
based  on  the  new  approach  by  the  finite  element  method,  for  an  example  of 
uniform  wall  thickness  of  the  nosed  shell.  It  was  demonstrated  that  both 
approximate  solutions  are  useful  in  design,  one  for  its  simplicity  and  the 
other  for  its  accuracy.  The  approximate  solutions  presented  here  will  also 
be  helpful  for  the  further  development  of  the  new  approach  to  preform  design 
in  shell  nosing. 


*  Part  2  of  "Preform  Design  in  Metal  Forming." 


Nomenclature 


r,  6,  z 


b,  B 


b 


£ 


» 


B 


£ 


£ 

t 

f 


-  radial,  circumferential  and  axial  coordinates  in  the  cylin¬ 
drical  system 

-  coordinate  along  the  die  contour 

-  angular  coordinate 

-  the  value  of  $  corresponding  to  the  nose  tip  where  the 
reaius  is  equal  to  b 

-  velocity  in  the  s-direction 

-  stress  components  in  the  s-direction  and  the  circumferential 
direction,  respectively 

-  flow  stress  (constant)  of  the  material 

-  strain-rate  components  in  the  s-direction,  the  circumferen¬ 
tial  direction  and  the  thickness  direction,  respectively 

-  radius  of  the  die  contour 

-  r,  z  -  coordinates  for  the  center  of  radius  of  the  die 
contour 

-  length  of  the  preform 

-  r  coordinate  and  corresponding  s  coordinate  at  the  entrance 
of  the  die 

-  tube  thicknesses  in  the  preform  and  in  the  nosed  shell, 
respectively 

-  radius  and  corresponding  s-value  at  the  nose  tip  of  the 
shell 

-  radius  and  corresponding  s-value  at  the  nose  tip  of  the 
shell  when  an  element  £  just  begins  to  deform 

-  axial  coordinate  in  the  preform 

% 

-  time 

-  subscript  referring  to  the  final  configuration 


1 .  INTRODUCTION 

Shell  nosing  is  the  process  of  forming  an  ogive  nose  at  an  open  end  of  a 
tubular  part  by  pressing  the  tube  into  a  suitably  formed  die.  The  internal 
cavity  of  a  tubular  part  (preform)  is  formed  to  finished  shape  through  a 
series  of  piercing  and  drawing  operations,  and  the  open  end  of  the  preform 
is  then  contoured  by  a  nosing  press,  as  shown  in  Fig.  1.  The  metal  flow  in 
nosing  is  complex.  In  general,  wall  thickness  increases,  but  the  shell  may 
elongate  or  shrink  in  length  after  the  operation,  depending  upon  friction. 
Understanding  the  deformation  mechanics  of  the  process  is  an  essential  step 
toward  the  control  of  metal  flow  and  proper  design  of  the  process. 

Because  shell  specifications  require  certain  wall  thickness  distributions 
after  nosing,  the  design  of  the  preform  with  which  the  desired  shape  after  the 
nosing  operation  can  be  achieved  is  a  very  important  problem.  For  certain 
shells,  uniform  wall  thickness  is  required,  whereas,  in  other  shells,  wall 
thickness  distribution  after  nosing  is  defined  by  two  ogive  radii. 

For  shells  with  uniform  wall  thickness  after  nosing,  Nadai  [1]  outlined 

an  approach  for  designing  the  preforms  based  on  the  strains  in  the  shell  nose. 

/ 

Using  Nadai' s  approach,  Carlson  [2]  suggested  the  method  for  determining  the 
original  shell  profile,  knowing  the  final  shape.  The  formulation  for  the  dis¬ 
tribution  of  thickness  in  the  preform  was  derived  from  the  geometry  of  de¬ 
formed  configurations  and  assumed  distributions  of  the  strains.  Thus,  the 
final  expression  for  the  determination  of  preform  thicknesses  requires  an 
estimate  of  a  quantity,  for  example,  elongation  at  the  tip  of  the  nose. 

Lahoti  et  al.  [3]  developed  the  computer  program  for  the  preform  design  pro¬ 
cedure  based  on  Carlson's  method.  Tang  et  al .  [4],  assuming  that  the  strain 
distributions  that  resulted  from  the  finite-element  simulation  are  valid  for 
preforms  having  nonuniform  wall  thickness,  designed  two  types  of  preform 
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shapes  to  obtain  uniform  wall  thicknesses  after  nosing.  In  this  paper  two 
approximate  solutions  for  designing  the  preforms  are  presented.  The  degree 
of  approximation  of  the  solution  is  examined  by  the  accurage  numerical  solu¬ 
tion  based  on  a  new  approach  by  the  finite  element  method  [5],  for  the  case  of 
uniform  wall  thickness  of  the  nosed  shell. 


2.  PRELIMINARIES 

Figure  2  shows  the  die  geometry  and  the  coordinate  system.  In  the  cylin 
drical  coordinate  system  (r,  z,  8)  with  the  origin  0  at  the  die  entrance,  the 
die  contour  is  expressed  as  r  =  R  si nd)-aQ ,  and  z  =  R  cos4>-bo.  The  coordinate 
s  measures  the  distance  along  the  die  contour  from  the  origin  O'  where  r  *  0. 
With  reference  to  Fig.  2,  and  using  the  notation  given  in  the  nomenclature. 


the  basic  equations  for  shell  deformation  are  given  by 
(i)  equilibrium  equation 


if  (h  o,»  ♦  7  <”s  -  °e>  ■  0 


(ii)  yield  condition  (plane  stress) 


2  A  2  2 
°s  ’  °s  °e  +  °e  =  °o 


(iii)  flow  rule 


2  °s  *  °e 


2  °8  '  °s  -(ce  *  “s' 


where  strain-rates  are  defined  as 
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(Vs  <  0) 


(1) 

(2) 

(3) 
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Since  the  thickness  h  in  general  is  a  function  of  coordinate  and 
time  t,  the  incompressibility  condition  is  written  as  [6,7] 


3 


I 


I 


1  Dh  , 
h  Dt 


3V  Vs 

—  cost!)  +  7-  cos<j) 


=  0 


“here  m 5  Ji  *  vs cos*  j? 


(5) 


Note  that  if  h  is  a  function  of  space  variables  only,  eq.  (5) 
is  reduced  to 

^  (h  r  Vs)  *  0  (5a) 

In  Fig.  3  the  initial  and  the  two  configurations  of  the  shell  nose  are 
shown.  Figure  3(b)  shows  the  configuration  of  shell  nose  at  a  time  when  the 
element  £  just  begins  to  deform,  and  characterized  by  the  radius  of  nose  tip 
denoted  by  b  =  b^.  Identifying  an  element  in  the  initial  and  the  deformed 
configurations,  the  incompressibility  equation  can  be  written  as 

h  <  r  •  ds  *  h.  •  rQ  •  d£  (b,  constant)  (5b) 

and 


h. 

h 


r  ds 

ro  . 


(6) 


When  the  final  configuration  is  specified,  the  radius  of  nose  tip  b  is  given 
and  h  is  known  as  a  function  of  r.  Thus  the  problem  essentially  is  to  determine 
the  relationship  between  the  corresponding  coordinates  in  the  initial  and 
final  configurations,  namely,  the  relations  between  5  and  s  (or  5  and  r). 


3.  APPROXIMATE  SOLUTION 
Procedures: 

Taking  the  magnitude  of  the  velocity  V$  at  the  nose  tip  as  unity,  the 
nose  tip  coordinate  B  (or  b)  can  be  used  as  the  time  scale  (dt  *  -  dB  * 
-db/cos4»h).  Then,  applying  the  condition  that  B  a  srt  at  E,  ~  0  and  B  ■  Br 

u  Os 


(7) 


(>(> 


at  £,  we  have 


<v 


dB 


s=s. 


From  the  definition  of  the  velocity  V  , 

V 

ivsi  ■  <i> 


(8) 


with  the  initial  condition  that  s  =  sq  at  B  =  B^.  The  subscript  £  in  eq.  (8) 
denotes  an  element  £.  Thus,  the  relationship  between  E,  and  s  (or  £  and  r)  can 
be  obtained  by  eliminating  B.  from  eqs.  (7)  and  (8),  when  the  velocity  V  is 

S  ^ 

known.  The  determination  of  the  velocity  V$  is  achieved  from 
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VC 

2o  -  o.  1  r  > 
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e  "s 

resulting  from  eqs.  (3)  and  (4). 

Integration  of  eq.  (9),  with  approximate  stress  solutions  [8],  given  by 

1,  results  in 


°s  r  °o 

=  -  Jin  £  ,  and  —  = 
o  b  CT 
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|*,|  *  <?>2  0  -  hn  S>'3 


(10) 


The  derivation  of  Eq.  (10)  is  given  in  the  Appendix.  In  order  to  perform  inte¬ 
gration  of  V$,  involved  in  Eqs.  (7)  and  (8),  further  approximation  of  the  velo¬ 
city  solution  (10)  is  necessary,  and  the  final  solution  for  pre-form  shapes  de¬ 
pends  on  this  approximation. 

Solution  I: 

A  simple  approximation  [9]  for  the  velocity  V$  is  given  by 

i*,i  *Vf  (in 

The  velocity  distributions  given  by  eqs.  (10)  and  (11)  are  compared  in  Fig.  4. 
Substituting  eq.  (11)  into  eq.  (7) 


■tm 


=  j=  W.)  -  «b5)] 

V  0 


G(x)  = 


M+an 


From  eq.  (8) 


.  dr  ,  JT 

db  cosd>  *  r 


Integrating  from  rQ  to  r  for  r  and  from  b^  to  b  for  b,  we  have 
G(r)  -  G(r0)  *  G(b)  -  G(b^) 

Eliminating  b^  from  eqs.  (12)  and  (13),  we  obtain 


[G(r)  -  G(b)] 


h,  dr  .  /  dG(r) 

From  eq.  (6)  T  •  —  ^-Qs0  for  b  *  constant,  but  ^  *  fT^/  dr 


V^o 

1  r* 


cos<b  from  eq.  (14). 


Thus  T’^TT  (,5) 

It  is  of  interest  to  note  in  eq.  (15)  that  h  is  not  a  function  of  time  and 
can  also  be  obtained  directly  from  the  incompressibility  condition  (5a)  with 
the  velocity  V$  given  by  eq.  (11).  Noting  that  usually  a0  +  r0  =  R  and 

cos<p  2  (1  -  ■■£-)  ,  the  expression  for  E  in  eq.  (14)  becomes 

s..frr  [.u-'Vf- ^fTT x-„.v-r,)frvrrL- 


The  derivation  of  eq.  (16)  is  given  in  the  Appendix. 
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The  solution  for  a  preform  shape  is  given  by  eqs.  (15)  and  (16).  The  quanti 
ties  h,  r  and  b  in  the  equations  refer  to  those  of  given  configuration  after 
nosing,  namely,  h^,  r^  and  b^  . 


Equation  (17)  is  also  compared  with  the  velocity  distribution  given  by  eq. 
(10)  in  Fig.  4.  It  can  be  seen  in  Fig.  4  that  this  approximation  is  better 
than  a  simple  one  used  for  solution  I.  In  order  to  perform  integrations  in¬ 
volved  in  the  solution  procedure,  further  approximation  must  be  made,  parti¬ 
cularly,  with  reference  to  the  relation  between  s  and  r.  It  is  assumed  that 


dr 

ds 


a  cos<|>  *  cosij»m  (=  const) 


(18) 


over  the  range  <  r  <  ro  where  4>m  is  an  average  value  of  q  over  the  rele¬ 


vant  range. 

With  this  assumption,  eq.  (7)  becomes 


C 


db 

cos<{>b 


If  we  define  an  average  value  ,  according  to 


For  derivation  of  eq.  (20),  see  the  Appendix.  Equation  (20)  gives  the  relation 
between  and  5  for  given  die  deminsions  and  nosed  shell  configuration. 
Equation  (8),  with  the  assumption  given  by  eq.  (18),  now  becomes 


<l>5  *  <&>e  *  i»,i  • 

and  its  integration  (given  in  the  Appendix)  results  in 


r  /T"  Jin  £  »  rQ  /"b^”  in  ^  (21 ) 

The  thickness  ratio  given  by  eq.  (6),  according  to  eqs.  (20)  and  (21),  becomes 


The  derivation  of  eq.  (22)  is  given  in  the  Appendix.  Equations  (20),  (21)  and 
(22),  if  we  use  hp  rf  and  b^  for  the  corresponding. quantity,  constitute  the 
solution  II  for  preform  design.  In  passing  it  is  to  be  noted  that  the  thickness 
change  expressed  by  eq.  (22)  resembles  that  of  the  flange  in  deep  drawing  obtain 
ed  by  Hill  [6],  It  should  be  further  noted  that  the  thickness  h  given  by  eq. 
(22)  is  not  only  dependent  on  r  but  also  a  function  of  tine  (b)  and  satisfies  in 
compressibility  condition  (5)  with  the  approximation  (18). 

4.  NUMERICAL  EXAMPLE  AND  DISCUSSION 

The  preform  dimensions  (h.  as  a  function  of  E)  are  calculated  using  the 
approximate  solutions  obtained  in  the  previous  section,  for  the  case  of  uni- 
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form  thickness  after  nosing. 

The  specifications  for  the  die  geometry  and  the  nosed  shell  dimensions 
(unit:  mm)  are  as  follows: 

R  =  259.00  (R/rQ  =  12.03)  (a  =  237.20,  bQ  =  18.75) 

rQ  *  21.53,  bf  *  14.47  (bf/rQ  =  0.671) 

hf  =  h0  *  6.90,  Lf  =  48.84  (Lf/rQ  =  2.27) 


The  computation  for  the  preform  dimensions  with  Solution  I  is  straight 
forward,  and  the  results  for  the  above  specifications  are  tabulated  in  Table 
1  and  shown  graphically  in  Fig.  5a. 


Table  1.  Preform  Dimensions  According  to  Solution  I. 


r/ro 

Vhf 

-(h./hf  -  1)  x  100 

0.671 

0.000 

0.819 

18.1 

0.700 

0.105 

0.837 

16.3 

0.725 

0.202 

0.851 

14.9 

.  0.750 

0.304 

0.866 

13.4 

0.775 

0.414 

0.880 

12.0 

0.800 

0.532 

0.894 

10.6 

0.825 

0.660  • 

0.908 

9.2 

0.850 

0.799 

0.922 

7.8 

0.875 

0.953 

0.935 

,  6.5 

0.900 

1.125 

0.949 

5.1 

0.925 

1.323 

0.962 

3.8 

0.950 

1.562 

0.975 

2.5 

0.975 

1.877 

0.987 

1.3 

1.000 

2.649 

1.000 

0 

The  preform  dimensions  according  to  Solution  II  are  calculated  by  determining 
bc/ro  for  given  r/rQ  from  eq.  (21)  first,  and  then  from  eqs.  (20)  and  (22). 
The  computed  results  are  given  in  Table  2. 


% 


Table  2.  Preform  Dimensions  According  to  Solution  II. 


r/ro 

Vro 

«/ro 

W 

-(hi/hf  -  1)  x  100 

■351 

0.819 

18.1 

mmm 

17.0 

HSU 

0.291 

0.841 

15.9 

0.431 

14.8 

0.574 

0.864 

13.6 

0.719 

0.877 

12.3 

HE IB 

0.866 

0.891 

10.9 

K'fVriSH 

1.016 

0.905 

9.5 

rliYfB 

1.167 

0.920 

8.0 

Kk  ft 

1.320 

0.935 

6.5 

■M 

1.475 

0.951 

4.9 

HUS 

0.728 

1.632 

0.967 

3.3 

0.700 

1.790 

0.984 

1.6 

1.000 

0.671 

1.950 

1.000 

0 

The  preform  dimensions  tabulated  in  Table  2  are  plotted  in  Fig.  5b. 

For  preform  design  in  shell  nosing,  the  accurate  solution  can  be  obtained 
by  a  new  approach  -  Backward  Tracing  Scheme  -  using  the  finite  element  method 
[5].  The  concept  involved  in  the  approach  is  to  trace  the  loading  path  back¬ 
ward  from  the  final  configuration  by  the  rigid-plastic  finite  element  method. 
The  two  preforms,  one  having  distributed  inner  radii  and  the  other  with  a 
cylindrical  inner  surface,  were  designed  for  the  present  example.  The  results 
are  almost  the  same  for  the  two  types  of  preforms,  and  the  preform  dimensions 
for  the  case  where  the  outer  diameter  of  the  preform  is  constant,  are  given 


in  Table  3, 
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Table  3.  Preform  Dimensions  -  Numerical  Solution. 


5/r°. 

Vf 

-(hi/hlf  -  1)  x  100 

0 

0.861 

13.9 

0.220 

0.846 

15.4 

0.435 

0.857 

14.3 

0.650 

0.878 

12.2 

0.863 

0.897 

10.3 

1.076 

0.917 

8.3 

1.288 

0.932 

6.8 

1.500 

0.952 

4.8 

1.710 

0.968 

3.2 

1.920 

0.974 

2.6 

2.130 

0.987 

1.3 

2.350 

0.996 

0.4 

The  degree  of  approximation  of  the  solutions  I  and  II,  therefore,  can  be 
examined  by  comparing  them  with  the  accurate  finite  element  solution.  A  com- 
parison  of  these  solutions  is  shown  in  Fig.  6.  It  is  seen  that  the  approxi¬ 
mate  solution  II  gave  an  excellent  estimate  of  preform  dimensions  for  the  most 
of  the  length,  except  near  the  preform  tip.  Approximate  solution  I  deviates 
from  the  numerical  solution  somewhat,  however,  simplicity  is  its  strength. 
Another  important  dimension  of  the  preform,  besides  the  thickness  distribution, 
is  the  total  length  of  the  preform.  In  Fig.  6,  (S/r0)0  is  the  length  of  the 
part  which  goes  into  the  die  after  nosing.  Thus  the  total  length  of  the  preform 
Lp  is  given  by  the  sume  of  (C)Q  and  the  constant  LQ. 

Solution  I  over-estimate  the  total  length,  but  solution  II  again,  is  very 
close  to  the  numerical  solution.  It  is  seen  that  in  the  numerical  solution 
the  thickness  at  (C/rQ)0  is  slightly  smaller  than  the  final  thickness.  This 
reflects  the  fact  that  the  actual  deformation  zone  extends  beyond  the  die 
entrance  during  nosing. 
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5.  CONCLUDING  REMARKS 

The  two  approximate  solutions  for  preform  den3n  in  shell  nosing  were 
derived,  and  it  was  shown  that  they  are  both  useful  in  design,  particularly, 
one  for  its  simplicity  and  the  other  for  its  accuracy.  The  examination  of 
approximation  was  possible,  because  of  the  new  numerical  approach  to  the  pro¬ 
blem  by  the  finite  element  method.  The  approximate  solutions  in  this  paper 
do  not  take  into  account  the  effects  of  interface  friction  and  material  pro¬ 
perties.  In  this  respect,  the  new  approach  will  assume  increasing  importance 
in  preform  design  in  the  future.  However,  the  approximate  solutions  developed 
here  will  be  useful  in  production  at  the  present  time.  Also,  it  is  antici¬ 
pated  that  the  approximate  solutions  can  play  a  role  during  the  course  of  further 
development  of  the  new  approach  as  an  initial  solution  from  which  the  final 
solution  will  converge  by  iteration. 
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APPENDIX 

Derivations  of  the  Equations 

( 1 )  Derivation  of  Eg.  (10). 

Substitution  of  approximate  stress  solutions  given  by 
°s  r  aA 

—  =  -  tn  r-  and  —  =  -  1 ,  into  equation  (9),  results  in 
°o  D  ao 

>  ’  (■  2  *  — - — p\  7- 

5  V  2  - ln  f/ 

Integration  can  be  performed  by  letting  tn  £  *  x  and  dx  *  ,  which 

gives 

tn  |  V  |  *  -  2  tn  £  -  3  tn  (2  -  tn  £■)  +  C 

From  the  boundary  condition  that  | Vg |  *  1  at  r  =  b,  the  integration  con¬ 
stant  becomes  C  ■  3  tn  2.  Then, 

tn  |V$|  *  -  2  tn  £  -  3  tn  (1  -  \  tn  , 

from  which  eq.  (10)  results; 

|VSI  -  (£)  (1  -  \  tn  £)'3 

(2)  Derivation  of  Eq.  (16). 

From  eq.  (14) 

- (G(r)  -  G(b)) 

0  ro  ^ 


where 


G(*>  .  f  —  - dx  =  f  gz  ■ 


Let  *  sin^y,  then  —  *  2  siny  cosy  dy,  and  G(x)  is  expressed  by 
ro  ro 


We  start  with 


where  cos<b  .  from  eq.  (19),  is  expressed  by 


15  77 


(4)  Derivation  of  Eg.  (21). 

In  order  to  integrate 


Since  the  subscript  £  denotes  an  element  £,  the  limits  of  integration 
are  b  *  when  r  *  rQ  and  b  =  b  when  r  *  r  for  the  element  £.  Thus, 


(5) 


which  gives  eq.  (21). 

Derivation  of  Eg.  (22). 

From  eqs.  (6)  and  (18),  we  have 

hi_  m  jr_  ds_  a  jr_  1  dr 

h  r  d£  r  cosa  cF  * 

o  ^  o  m  ^ 


for  b  »  const. 


» 


where 

dr  .  dr  dbS 

I  if  • 

Differentiating  eq.  (21)  with  b  *  constant,  we  obtain 


dr  _  ro 


FIGURE  TITLES 


Figure  1.  Shell  Nosing  Process  (from  Ref.  [3]). 

Figure  2.  Die  and  Workpiece  geometries  and  the  coordinate  system. 

Figure  3.  3  Stages  of  the  Shell  configurations. 

Figure  4.  Comparison  of  various  velocity  distributions. 

Figure  5.  Preform  dimensions,  based  on  (a)  solution  I  and  (b)  solution 
II. 


Figure  6. 


Comparison  of  the  approximate  solutions  with  the  solution 
by  the  new  approach  using  the  finite  element  method. 
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Comparison  of  the  approximate  solutions  with  the  solution 
by  the  new  approach  using  the  finite  element  method. 


